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ABSTRACT 


This  report  is  primarily  concerned  with  the  stud>  of  t 

electromagnetic  scattering  by  random  scatterers.  Potential 

s 

applications  to  radar  detection  and  communication  problems 
are  stressed.  In  radar  detection  problems,  it  is  often  necessary 
to  detect  a  target  echo  in  the  presence  of  other  unwanted  echoes 
(clutter).  In  order  that  a  radar  receiver  can  be  designed  to  operate 
effectively  in  the  presence  of  clutter  interference,  it  is  necessary 
to  develop  a  suitable  theoretical  model  of  the  clutter  which  can  be 
used  for  the  design  and  evaluation  of  detection  schemes.  In  radio 

J 

communication  problems,  it  is  desirable  that  signals  can  be 

*  .. 

communicated  between  two  non-line-of- sight  points  by  means  of 
electromagnetic  scattering  from  a  medium  which  occupies  the 
region  of  space  illumi'  ited  by  the  two  antenna  beams.  Detailed 
knowledge  of  the  scattering  characteristics  of  the  medium  will 
enable  one  to  select  proper  signal  parameters  for  transmittal  of 
information  and  optimum  processing  schemes. 

A  cloud  of  dipoles  (chaff),  dispensed  in  a  proper  region 
in  space  to  act  as  radar  reflectors,  can  be  described  as  an 

* 

assembly  of  random  scatterers.  When  chaff  dipoles  are  dispensed 

from  a  moving  craft  in  space,  they  will  in  general  move  relative  * 

ii 


to  one  another.  In  addition,  each  dipole  will  have  a  tumbling 
motion  due  to  effects  of  injection  forces,  body  instability,  and 
other  aerodynamical  properties.  Since  the  location  and  motion 
of  the  individual  scatterers  are  unknown,  the  scattering  problems 
are  best  treated  statistically. 

The  bistatic  radar  reflecting  characteristics  of  a 
randomly  tumbling  dipole  are  investigated.  An  expression  for 
the  scattered  voltage  is  derived  by  application  of  the  Lorentz 
reciprocity  theorem.  The  correlation  properties  of  the  received 
signal  are  examined.  Some  statistical  assumptions  are  made  in 
order  to  obtain  readily  usable  results. 

A  theoretical  model  is  developed  for  the  radar  echo  from 
a  random  collection  of  moving  dipole  scatterers.  The  analysis 
of  the  model  takes  into  account  some  effects  of  scatterer  rotation 
which  have  been  neglected  in  previous  work.  The  fluctuating 
characteristics  of  clutter  echoes  are  also  determined.  The 
theory  and  some  experimental  results  in  the  literature  are  shown 
to  be  in  relatively  good  agreement. 

The  properties  of  random  scatter  communication  channels 
are  also  investigated.  The  constitutive  parameters  of  the 
scattering  medium  are  assumed  to  be  varying  randomly  with 
space  and  time.  The  effects  of  antenna  gain  are  included  in  the 
derivation  of  the  channel  function  in  order  to  take  explicit  account 
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Chapter  1 

4  INTRODUCTION 

1.  1  The  Problem  and  Objectivea  of  thia  Study 

In  radar  system*  it  is  often  necessary  to  detect  a  target  echo 
in  the  presence  of  other  unwanted  echoes  or  clutter.  Clutter  echoes 
include  signals  reflected  from  chaff,  surface  of  the  ground  end  sea* 
vegetation,  etc.  The  ability  of  a  radar  to  detect  a  target  echo  is 
limited  by  the  presence  of  clutter.  In  order  to  develop  a  radar 
receiver  that  will  operate  effectively  in  the  presence  of  clutter  later-' 
ference*  it  is  desirable  to  have  a  theoretical  model  of  the  clutter 
1  which  can  be  used  for  the  design  and  evaluation  of  detection  schemes. 

If  the  clatter  waveform  were  precisely  known,  it  would  be  possible* 
in  principle,  to  synthesize  a  filter  which  optimizes  the  signal-to- 
noise  ratio  at  the  output  o'  the  receiver. 

The  clutter  model  most  often  assumed  for  theoretical  aoalys#s 
is  the  "point  scatterer"  model.  Rotational  motion  of  the  scatterers 
is  neglected.  The  scatterers  are  treated  as  point  targets  with 
variable  cross  section. 
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Chaff  is  a  form  of  countermeasure  used  against  radar.  A 
chaff  cloud  usually  consists  of  a  large  number  of  dipoles  resembling 
an  assembly  of  random  scatterers.  When  chaff  dipoles  are  dispensed 
from  a  moving  craft  in  space,  the  dipoles  experience  certain  forces 
which  may  cause  them  to  rotate  or  tumble.  In  some  applications  As 
dipoles  may  be  constrained  tc  rotate  "bout  a  preferred  axis.  Thus, 
ignoring  the  effects  of  scattefer  rotation  will  result  in  an  incorrect  ' 
representation  c  the  model.  Ground  clutter  exhibits  echoing  charac- 
teristics  similar  to  those  of  an  assembly  of  random  scatterers. 
Rotational  motion  of  the  scatterers,  or  equivalent  movement®  of 
branches,  leaves,  grass,  etc. ,  under  the  effects  of  wind  force,  can 
become  a  major  contributing  factor  to  the  observed  fluctuations  of 
echo  intensity. 

Herein,  a  new  theoretical  model  for  the  radar  echo  from  such 
collections  of  random  scatterers  is  developed.  The  scatterers  are 
treated  as  randomly  oriented  dipoles  instead  of  points  with  variable 
cross  sections.  This  assumption  permits  the  determination  of 
scatterer  rotation. 

In  many  radio  communication  problems  signals  must  bo 
communicated  between  two  non-line-of-sight  points  by  means  of 
electromagnetic  scattering  from  a  medium  which  occupies  the  region 
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of  spice  illuminated  by  the  two  antenna  beams.  The  scattering 
medium  may  be  continuous  such  as  the  inhomogeneous  atmosphere  or 
may  consist  of  a  random  collection  of  metallic  scatterers  such  as 
chaff. 

The  scattering  medium  in  general  varies  with  both  position 
and  time.  In  the  analyses  of  random  scattering  most  investigators 
have  ignored  the  time  dependence  and  have  been  concerned  only  with 
the  calculation  of  the  average  cross  section  of  the  scattering  medium. 
The  average  cross  section  is  useful  for  estimating  the  signal-to -noise 
ratio  at  the  receiver;  however,  it  does  not  provide  adequate  informa¬ 
tion  regarding  the  fluctuating  characteristics  of  the  scattered  wave- 
form.  In  the  present  work,  we  study  scatter  channels  in  which  the 
scattering  medium  is  randomly  varying.  General  time  dependence  is 
assumed.  The  statistical  nature  of  such  channels  are  carefully 
examined. 

1.  2  Summary  of  Previous  Work 

Clutter  echoes  are  often  regarded  as  signals  reflected  from 
an  assembly  of  independent  random  scatterers.  A. cloud  of  dipoles 
(chaff),  dispensed  randomly  in  space,  is  an  example  of  such  a  set. 
Because  of  the  r;tndom  location  and  motion  of  the  individual  scatter¬ 
ers,  clutter  echoes  arc  best  treated  statistically.  Siegertf26], 
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Kerr  [13],  tod  Lawton  and  Uhlenbeck  [27],  developed  the  first 
probability  distribution  functions  for  the  magnitude  of  the  returned 
ecliJ  from  an  infinitely  dense  chaff  cloud.  Kelly  and  Lerner  [24] 
developed  a  more  extensive  theoretical  chaff  model.  They  showed 
that  the  returned  echo  from  the  chaff  cloud  could  be  regarded  as  & 
random  process  and  treated  the  scatterers  as  points  with  variable 
cross  sections, 

Kerr  [13]  and  Lawson  and  Uhlenbeck  [27]  presented  some 
experimental  results  on  the  measurements  of  chaff  models  and  ground 
clutter.  Some  anomalous  characteristics  were  observed  in  the 
measured  correlation  functions.  Kerr  suggested  that  the  cause  could 
have  been  due  to  the  effects  of  scatterer  rotation,  but  he  made  no 
theoretical  analysis  to  clarify  this  phenomenon. 

The  possibility  of  using  a  cloud  of  chaff  dipoles  for  estab- 
fishing  short  time  communications  between  non-line -of- sight  points 
was  considered  by  Blom  [8]  and  Hessemer  [9].  The  recent  West  \ 
Ford  Project  [11]  extended  this  basic  c  oncept  and  deployed  an 
orbiting  dipole  belt.  This  belt  established  the  first  world-wide 
scatter  communication  channel. 

Some  of  the  statistical  properties  of  random  scatter  channels 
were  presented  recently  by  Kelly  [14].  Kelly's  analyses  include  both 


the  continuous  And  discrete  scatterer  channels.  He  took  into  account 
quite  generally  the  varying  character  of  the  constitutive  parameters 
of  the  medium  with  position  and  time. 

1.  3  Summary  of  Results 

Chapter  2  presents  an  analysis  of  the  reflecting  characterise 
tics  of  a  rotating  dipole.  The  dipoie  is  assumed  to  be  rotating  about 
an  arbitrary  axis  perpendicular  to  its  length.  Arbitrary  polarisation 
is  assumed  for  both  the  transmitting  and  receiving  antennas.  The 
Lorentz  reciprocity  theorem  is  applied  to  derive  an  expression  for 
the  bistatic  scattered  voltage.  It  is  shown  that  the  received  signal  is* 
in  general  modulated  in  both  amplitude  and  phase.  The  correlation 
properties  of  the  received  signal  are  examined.  Some  statistical 
assumptions  are  made  in  order  to  obtain  tractable  results. 

•  i 

A  theoretical  model  for  the  radar  echo  from  a  random  collec- 
ticm  of  dipole  scatterers  is  fomulated  in  Chapter  3.  The  analysis 
of  \the  model  takes  into  account  the  effects  of  scatterer  rotation.  An 
expression  for  the  correlation  function  of  the  returned  echo  is  derived 
in  terms  of  the  characteristics  of  the  transmitted  waveform,  polari¬ 
zation,  and  the  statistical  properties  of  the  scatterers.  Under  ths 
assumption  that  the  cloud  density  is  slowly  varying,  the  process 
becomes  stationary.  The  probability  density  of  the  returned 


waveform  ii  shown  to  approach  a  Gaussian  probability  density  as  the 
rate  of  echo  retorn  becomes  'large.  The  theory  can  be  applied  to  the 
study  of  a  wide  class  of  clutter  signals  which  exhibit  similar  echoing 
characteristics.  For  example,  targets  making  up  ground  clutter  such 
as  trees,  grass,  etc. ,  will  move  and/or  oscillate  with  the  wind.  The 
individual  blades  or  stems  scatter  like  dipoles.  The  fluctuating 
characteristics  of  clutter  echoes  are  also  determined.  Comparison 
between  the  theory  and  some  experimental  results  in  the  literature 
are  shown  to  be  in  relatively  good  agreement. 

Chapter  4  describes  a  communication  channel  which  is 
provided  by  means  of  electromagnetic  scattering  from  a  collection  of 
scatter*.3*  randomly  distributed  in  a  region  of  space.  The  collection 
of  scatterers  behaves  like  a  dispersive  medium,  and  the  number 
density  varies  with  boih  position  and  time.  The  channel  response 
function  is  described  in  terms  of  the  scattering  properties  of  the 
medium.  The  mean  and  the  covariance  of  the  channel  output  are 
derived  under  the  assumption  of  Poisson  statistics.  The  gain  func¬ 
tions  of  the  transmitting  and  receiving  antennas  are  included  in  the 
analysis  in  order  to  take  explicit  account  of  the  iact  that  scatterers 
may  flow  in  and  out  of  the  volume  illuminated  by  the  two  antenna 
beams.  The  effects  of  acatterer  rotation  are  also  considered.  The 

v 

case  of  dipole  scattering  is  given  as  an  example. 


The  properties  of  a  dispersive  continuous  channel  are 
considered  in  Chapter  5.  The  scattered  fields  are  derived  in  terms 
of  the  electric  and  magnetic  susceptibility  functions  of  medium.  As 
in  the  discrete  scatterer  case,  the  effect  of  antenna  gain  is  also  taken 
into  account  in  deriving  the  channel  function.  Since  there  are  many 
similarities  in  the  formalism  between  the  two  types  of  scattering, 
most  of  the  results  obtained  for  the  discrete  scatterer  case  are  valid 


for  the  continuous  case.  * 


Chapter  2 


BISTATIC  RADAR  REFLECTIONS  FROM  A 
RANDOMLY  ROTATING  DIPOLE 

2. 1  Introduction 

'When  a  plane  electromagnetic  wave  impinges  upon  a  target 
in  space,  some  of  the  incident  energy  will  be  scattered  in  the  direc¬ 
tion  of  the  radar  receiver.  The  characteristics  of  the  scattered  • 
signal  depend  upon  (1)  the  size,  shape,  and  orientation  of  the  target 
and  (2)  the  frequency  and  state  of  polarization  of  the  incident  wave. 

The  problem  of  electromagnetic  scattering  from  cylindrical 
wires  has  been  considered  by  a  number  of  investigators  [1]  -  [7]. 

A  bundle  of  metallic  wires  or  strips  (known  as  chaff),  cut  to  resonate 
at  a  certain  wavelength,  may  be  dispensed  in  a  proper  region  in 
space  to  act  as  effective  radar  reflectors.  These  reflectors  can 
perform  two  different  functions: 

(1)  To  present  certain  characteristic  echoes  to  a  radar  for 
purposes  of  identification,  deception,  or  confusion  [3]. 

(2)  To  provide  a  channel  for  beyond  line-of-sight  communi¬ 
cations  [8]  •  [11]. 
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A  cloud  of  dipolus  can  be  represented  by  an  assembly  of 
random  scatterers.  When  chaff  dipoles  are  dropped  info  space,  they 
will  in  general  move  relative  to  one  another.  Each  dipole  will  have 
a  slightly  varying  position  and  velocity,  thus  giving  rise  to  a  Doppler 
spread  in  the  scattered  signal.  In  addition,  the  dipoles  will  have  a 
tumbling  motion  due  to  effects  of  body  instability  and  other  aero* 
dynamical  properties.  Tumbling  motion  of  the  dipoles  can  cause 
further  frequency  spreading  in  the  received  waveform.  In  many 
practical  situations,  tumbling  speeds  were  not  thought  to  be  very  * 
rapid  [11]  *[12]  and  considered  to  give  only  second  order  effects; 
while  in  other  situations,  the  effects  of  scatterer  tumbling  (or  - 
rotation)  could  be  observed  and  found  to  have  significant  contributions 
upon  the  received  waveform  [1?].  It  will  be  shown  in  Chapter  1' 
that,  in  general,  the  effects  due  to  scatterer  rotation  cannot  be 
neglected  if  the  Doppier  spread  arising  from  the  velocity  fluctuations 
of  the  scatterers  is  small. 

In  order  that  a  radar  receiver  can  be  designed  to  operate  Zm 
an  optimum  fashion,  an  adequate  knowledge  of  the  amplitude  sad 
phase  characteristics  of  the  scattered  waveform  will  be  required. 

In  the  analyses  of  random  scattering,  most  investigators  ignore  the 
effects  of  scatterer  rotation  and  assume  a  "point  scatterer"  model, 
using  the  information  on  the  average  cross  section  of  the  individual 


io 


scatterers.  The  "point  scatterer"  model  is  valid  for  rotationally 

symmetric  scattcrers.  However,  such  a  representation  Is  not  ‘ 

*d»qur.te  for  scattered*  of  other  shapes  such  as  dipolos. 

This  chapte;  presents  an  analysis  of  the  reflecting  charac¬ 
teristics  of  a  tumbling  dipole.  It  is  assumed  that  the  tumbling. motion 
of  the  dipole  can  be  described  by  a  rotation  about  an  axis  perpendicu¬ 
lar  to  its  length.  An  expression  for  the  bistatic  scattered  voltage  is  - 
derived  as  a  function  of  the  transmitter  and  receiver  polarisations 
and  the  orientation  of  the  rotation  axis  of  the  dipole.  The  statistical 
properties  of  the  signal  reflected  from  a  collection  of  randomly 
rotating  chaff  dipoles  are  discussed  by  means  of  a  mathematical 

»  ► 

model.  The  correlation  function  of  the  received  signal  envelops  is 
derived  in  terms  of  the  charccteristics  of  the  transmitted  waveform, 
polarisation,  and  the  statistical  properties  of  the  scatterers. 

2.  2  Characteristics  of  ihe  Scattered  Signal 

An  expression  for  the  scattered  signal  can  be  derived  by 
application  of  the  Lorentz  reciprocity  theorem  (Appendix  1).  Tbs 
radar  system  consists  of  a  transmitter  which  illuminates  the  target 
and  a  receiver  at  which  we  evaluate  the  scattered  signal.  The  total 
voltage  received  is  the  sum  of  the  voltages  induced  by  the  transmitter 
field  and  the  field  scattered  by  the  target.  In  the  following  derivation. 


11 

we  shall  ignore  the  coupling  effects  between  the  transmitter  and 
receiver,  but  allow  for  variable  polarisations.  The  received  voltage 
is  then  given  by 

*  *  mT  J  V  l4v  <*•*» 

R  V 

•  • 

where  I_  is  the  feed-point  current  which  would  have  to  be  supplied  to 
the  receiving  antenna,  if  it  were  used  as  a  transmitting  antenna,  to 
create  the  field  E  in  the  absence  of  the  target,  and  J  is  the  current 
density  associated  with  the  dipole.  If  the  dipole  is  thin  and  If  the 
field  can  be  considered  uniform  throughout  the  small  volume 
containing  the  dipole,  the  integral  on  the  right  side  of  (2. 1)  can  be 

expressed  in  terms  of  the  induced  moment  on  the  d’ipole.  Thus, 

•  • 

J  Er*  J  dV  =  taJJ  *Er  (2. 2) 

V 

s 

The  induced  dipole  moment  can  be  written  as 

p  =  *  6  (d  .  E_)  d  (2.  3) 

r  ©  e  .  T 

where  E_  is  the  transmitted  electric  field,  <X  is  the  polarisability 
T  6 

1  * 

of  the  dipole  ,  and  d  is  a  unit  vector  parallel  to  the  dipole  axis.  The 

*Here  the  polarizability  is  treated  as  a  seals  ,  more  generally,  It 
is  a  tensor  CM]. 
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field  quantities  E_  and  E_  can  be  written  as  E_  «  E__e_  and 
AH  T  OX  X 

8  where  e^  and  eR  are  the  polarization  vectors  of  the 

transmitter  and  receiver*  respectively*  and  the  symbol  -  is  used  to 
indicate  a  complex  quantity.  Equation  (2. 1)  becomes  . 


is  t  8  <w 

o  o  e  *  .  *  *  • 

R  EOTEOR«d'*T»<^*R» 


2  1 


(2.4| 


If  the  same  antenna  is  used  for  transmitting  and  receiving 

m  m  »  **»  **» 

(monostatic  case),  we  have  eT»  eR*  e  and  EQT*  EQR*  (2.4) 

reduces  to 


is  (8  -  _ 

o  o  e  2  2.?  -.2 

.  *  -jj—  eo  <a..» 


(2.5) 


The  bistatic  scattering  geometry  is  shown  in  Fig.  1.  The 
transmitter  and  receiver  Are  separated  by  the  angle  0.  The  incident 
wave  is  assumed  to  be  propagating  along  the  positive  z-axis  of  the 
transmitter  (unprimed)  coordinate  uystem*  and  the  scattered  wave*  is 
observed  in  the  direction  of  the  negative  z'-axis  of  the  receiver 
(primed)  coordinate  system.  In  order  to  compute  the  bistatic 
scattered  voltage*  it  is  necessary  to  express  the  dipole  orientation 
and  the  polarization  vectors  in  one  common  coordinate  system.  For 
simplicity,  we  choose  the  bistatic  angle  0  such  that  it  represents  a 
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rotation  in  the  xz -plane.  Hence,  the  unit  vectors  (s' ,  t',  a* )  in  the 

®  y  ® 

receiver  coordinate  system  are  related  to  the  unit  vectors 

AAA 

(a  ,  a  ,  a  ) in  the  transmitter  coordinate  system  by  the  rotational 
n  y  _  a 

transformation. 


cos  3 

0 

sin  3  \ 

0 

1 

• 

(2.6) 

-  sin  3 

0 

cos  3  / 

Figure  2  shows  the  coordinate  system  for  the  rotating  dipole. 
For  the  present  purpose,  the  distance  between  the  dipole  and  the 
radar  is  assumed  to  be  fixed.  The  dipole  is  rotating  about  an  arbi¬ 
trary  axis  s  defined  by  the  angles  *  and  T).  Let  (u,  v,  s)  be -a  set  of 
orthogonal  axes,  then  u  and  v  lie  in  the  plane  of  rotation  of  the  dipole 
If  'J>r  is  the  angular  rotation  frequency  of  the  dipole  and  if  a  is  its 
initial  position  with  respect  to  the  u-axis,  then  the  instantaneous 

A 

position  of  the  unit  vector  d  may  be-  represented  by 


d  *  a  cos  ♦  +  a  sin  f 
u  v 


<2.  V 


A  A 

where  Q  =  uu  t  +  a  and  a  and  a  are  unit  vectors  along  the  u  and  v 
r  «  •  v 

axes,  respectively.  For  convenience,  the  u-axis  is  chosen  to  bo 
coincident  with  the  intersection  of  the  xy-plane  and  the  plane  of 
rotation  of  the  dipole;  then,  by  means  of  a  transformation  between 
the  (x,  y,  z)  and  (u,  v.  s)  coordinate  systems  [IS],  we  obtain 


d  s  L -  *  (»iR  H  coal  ♦  coa  5  cos  »|  slat) 

X 

*  *Jco*  1  cost  -  cos  5  Rio  T)  slot  I 

T 

♦  s#(*ia  §  si**)]  f2.«) 

•  A  A 

where  (t  ,  '  a  )  are  unit  vectors  along  the  (x,  y,  s)~&xe»,  resptc- 
lively. 

The  polarization  vector  of  an  arbitrary  polarized  plane  wave 
can  be  written  ss  the  sum  of  a  right-handed  circularly  polarized  wave 
and  a  left-handed  circularly  polarized  wave.  Then,  for  the  trans¬ 
mitter  polarization,  we  writ* 


and  a^  and  a^  are  complex  constants  satisfying  the  relation, 

is, i*  ♦  is/  •  i 

Thus  the  state  of  polarization  of  the  transmitted  plane  wave  is  com¬ 
pletely  determined  once  the  complex  quantities  and  a^  are  known; 


e,  g. ,  for  circular  polarisation,  we  have  either  (a^  *  1,  a^  *  0)  or 

***  m*  m* 

(a^  =  0,  a^=  l)  and  for  linear  polarisation,  we  have  a^*  a|  • 

^o 

( !  //2 )  e  where  <p  is  the  angle  between  the  electric  field  vector  and 

o 

the  x-axis  as  shown  in  Fig.  1.  Similarly,  the  receiver  polarisation 


can  be  expressed  in  terms  of  the  primed  coordinate  system.  Time.. 


1’ 

T  1 


(2.1fl> 


where 


~(a*  -  is*) 
/TV*  T 


k  *  —  (a*  ♦  ii’) 

n>- *  ' 


l  I*  * « 


Substituting  qi  (2.C),  (2.9),  and  (2. 10)  into  (2.4)  together  with  the  aid 
of  (2.6)  yields,  after  some  algebraic  manipulations. 


a  *  a 


[C  +  U 


r,  -«#■ 


e  'He 


(2.111 


where 


to  <  a 

0^  A  M  a  AS  SNS 

«. *  -7^  eoteor 


(2*  12} 


stwws»«taf 
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C  *  -  J  (  A  A'  -  BB'|  (2.  13) 

***  1  ^  a  ’ 

U  s  *  4  *  A  4  B){ A  ♦  B  )  (2.  14)  / 

t.  .ia.iiix'.vt  (2.isi 

A  »  —  C?..*1'1- (2.UJ 

n 

*  1  ,a»  +  i  n  <«  -in, 

B=  - U.e  *+ae  *)  cos  5  (2.175 

a  *  * 

+  O'* *‘“Z  I  '  *r co,J f )  ** * 11  ]  ****** 


*'  *  ~/l  {[(*1C°‘2  2  '  *r*i‘,22)*M'' 

-  *arCO,^)e J  coe  ^ 

*  [(*i  +  *r)  ®in  P]  ,in  (2*1*) 

For  the  case  of  monostatic  radar  (8  =  0),  i.  e. ,  the  same  antenna  is 
used  for  transmitting  and  receiving,  we  have  E^»  ®*OTf 
and-B*  =  B»  The  quantity  a^  may  be  interpreted  as  the  maximum 
scattered  voltage  with  both  the  transmitter  and  receiver  polarisation 
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vectors  parallel  to  the  dipole  axis* 

The  first  term  in  (2. 1 1)  is  independent  of  f  and  hence  its 
amplitude  does  not  vary  with  time.  The  received  signal  consist*  of 
three  spectral  components,  namely  a  carrier  component  plus  the 
upper  and  lower  sideband  components  separated  from  the  carrier  by 
twice  the  rotational  frequency  of  the  dipole.  Furthermore,  the 
amplitudes  of  the  upper  and  lower  sideband  components  are  in  general 
not  equal,  depending  on  the  radar  polarization  and  the  orientation  of 
the  rotation  axis.  Thus,  both  amplitude  and  phase  modulation  of  the 
scattered  signal  ean  result. 

Equation  (2.  11)  is  a  general  expression  for  the  bistatic  scat- 
tered  voltage  with  arbitrary  transmitter  and  receiver  polarizations. 
As  a  special  example,  consider  the  case  of  backscattering  of  a 
circularly  polarized  wave.  Assume  that  both  the  transmitting  and 
receiving  antennas  are  polarized  in  the  same  sense  and  that  ths 
rotation  axis  is  the  positive  z-axis  (clockwise  rotation  about  the  line 
of  sight).  Using  Eqs.  (2. 11)  to  (2. 19)  and  setting  t  *  uJ^t  +#,  ws 
find 

;  _  -«(.  «♦«»  . 

for  the  right-handed  circular  polarization  (a  =  a'  *  1  and  a,  *  a'  *  0), 

,  r  r  •  ® 

and 
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“  *iV  «*.**> 

for  left-handed  circular  polarisation  (a.  *  a*  *  1  and  a  *  ^  «  0L 

*  *  r  r 

Since  the  scattered  voltage  is  a  real  function  of  time,  we  multiply 

♦  is  ft 

(2.20)  and  (2. 21)  by  e  °  and  evaluate  the  real  part.-  Menca,  vs 
obtain' 

/  J  m. 

12  I*J  co*  tl*0  ■  2»r>  *♦  YR3 - -RHCP 

j  |aj  cos[(ao+2»^)t4YL]— LHCP 

•  • 

where  YR  and  Yj^  are  constants,  is  the  angular  carrier  frequency 
of  the  transmitted  wave,  and  x^  is  the  rotation  frequency  of  the 
dipole.  Thus,  the  original  transmitter  frequency  is  shifted.  The 
frequency  shift  is  either  downvard  or  upward  depending  on  whether 
the  polarization  is  right-handed  or  left-handed,  respectively.  More 
precisely,  the  frequency  of  the  returned  signal  is  shifted  down  if  the 
polarization  vector  and  the  dipole  are  rotating  in  the  sam«  dire  "ns 
and  is  shifted  up  if  they  are  rotating  in  opposite  directions.  This 
result  has  also  been  shown  by  Mikulski  [  16]  and  Montgomery  t  JTj 
by  means  of  periodic  scattering  matrices. 


i  i 


-t'  tO**^*'- 


v^»-rH  r;  "t-  .7.  /!* 
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2.  3  The  Bistatic  Scattering  Cross  Section 

An  expression  for  the  bistatic  cross  section  can  be  derive-*  bp 
considering  the  power  flow  in  the  transmitted  and  the  received  waves* 
Thus,  the  average  power  received  is 


*  dr  rf 

o 


(2.22) 


where  Z  is  the  characteristic  impedance  of  the  transmission  line, 
o 

Substituting  (2.4)  into  (2.  22),  we  obtain 


P  = 


"°  V  |EOTl*|E„R|2|(d..T)(d..,,)|* 
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OR' 


(2.23) 


SI 


Let  P^,  be  the  average  power  supplied  to  the  transmitting  antenna  to 
• 

create  the  field  i*i  free  space  and  let  be  the  gain  function* 

then  the  time  average  magnitude  of  the  Poynting  vector  at  a  distance 

r  from  the  transmitter  is 
T 


iSTl 


G  P 
T  T 

4 n  r. 


=  -  sPr^r 

2  M  '  OT* 


(2.24) 


Similar  definitions  can  be  applied  to  the  receiving  antenna.  Thus, 


.r  .  GRPR  1  ~\/o  is  i2 

•SR>  S  _ T  *2  V  5“  IEo*» 

4n  r_  o 


(2. 2$) 


wh«r«  PR  «  —  is  the  power  which  must  be  supplied  to  the 

receiving  antenna  in  order  to  create  the  field  15_De_  in  free  space. 

OR  R 

Combining  (2.  23), .  (2.  24),  and  (2. 25),  we  obtain 


P  = 
s 


2  2 
n  a 

e 


G  G_  P_ 
T  R  T 

ia*\Z  2  2 

<4">  rT  ,R 


|(d.ST)(3.eR)  y 


(2.24) 


If  Eq.  (2.  26)  is  compared  with  the  bistatic  radar  equation. 


we  find. 


P 

s 


°t°rV* 

2  2 

(4n)  rT  rR 


o 


(2.27) 


0  =  «  ko*a.2  M  'rni.ij?  (2.2») 


where  a  is  the  bistatic  cross  section  and  k^  =  2n/\  is  the  free  space 

wave  number.  The  quantity  k*<X  ^  /4n  is  just  the  broadside  cross 

o  e 

A  * 

section  of  the  dipole  when  the  incident  electric  field  vector  is  parallel 
to  the  dipole  axis  [l].  For  a  thin  dipole,  the  polarizability  is 
approximately  given  by  [18], 


where  t  and  a  are  the  length  and  radius  of  the  dipole,  respectively. 
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This  approximation  is  valid  when  »  1.  Using  (2.  II),  (2.  28)  may 
be  written  as 


Equation  (2.  29)  also  applies  to  small  resonant  (tuned)  dipoles 

and  is  approximately  valid  for  half-wave  resonant  dipoles,  since  As 

difference  between  the  patterns  of  a  half-wave  dipole  and  a  short 

dipole  is  small.  However,  the  quantity  0^  must  be  replaced  appro* 

•  * 

priately;  for  example,  the  broadside  cross  section  of  a  small 
resonant  dipole  is  approximately  given  by  [18), 

* 

aQ  ~  X2  *  .716  X1  (2.31) 

and  that  of  a  half-wave  resonant  dipole  is  [5]  -  [6]» 

a  mt  . 86X2  (2.32) 

o 

Of  particular  interest  is  the  computation  of  the  average 
bi static  cross  section  of  a  randomly  rotating  dipole.  Assume  that 
the  initial  orientation  of  the  dipole  axis  in  space  is  uniformly 
distributed  and  that  every  orientation  of  the  rotation  axis  in  space  is 
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equally  probable,  then,  from  Eq.  (2.  29),  we  have 


o  2n  n 


<0>  *  4?  J  J  <C2+  U2*  L2)  sin  ?  d?  dq 


o  o 


*  0  <C24U24L2> 

O 


(2.  33) 


where  C,  U  and  L  are  the  magnitudes  of  2,  0,  andT>,  respectively, 

and  the  symbol  <  >  is  used  to  denote  an  ensemble  average.  If  we 

2  2  2  2 

denote  the  quantity  <  C  +  U  +  L  »  by  <  b  >  .  (2.  33)  can  bo  written 

© 

as 


< 0>  =  o  < b  > 
o  o 


(2.  34) 


Using  Equations  (2. 13)  through  (2. 19),  we  find 


<b2>  = 
o 


(2.35) 


1  i  .  ~  ~  ^  ^  ^  ,7  V 

+  2  l^ai"*r^*i"  a'r>  *  tai  +  ar^*i  +  ar^  coa  >31*} 


Equation  (2.  35)  is  a  general  expression  for  the  average  rela¬ 
tive  bistatic  cross  section  with  arbitrary  transmitter  and  receiver 
polarizations.  Hence,  if  the  radar  system  is  equipped  with  polari¬ 
zation  diversity  capability,  the  polarization  of  the  transmitter, 
receiver,  or  both  can  be  varied  to  achieve  either  maximum  or 


V 


.  3ti. 
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minimum  cross  section  as  desired.  For  the  case  of  backscattering, 
the  average  cross  section  becomes 


(2. 36) 


Since  the  process  is  random,  the  measured  cross  section 
usually  fluctuates  about  the  average  value.  The  probability  density 
functions  for  the  radar  cross  section  of  N  randomly-oriented  dipoles 
have  been  discussed  by  Borison  [19],  Fielding  [20],  and 
Rheinstein  [21], 


2. 4  Signal  Correlation  Properties 

■The  average  cross  section  provides  useful  information  for 
computing  the  expected  signal  to  noise  ratio  at  the  receiver  [11], 
[22].  However,  knowledge  of  the  average  cross  section  alone  is  not 
sufficient  for  optimum  radar  receiver  design  since  only  very  little 
information  is  provided  regarding  the  fluctuating  character! sties* of 
the  received  waveform.  Proper  design  of  radar  receivers  re^elres 
a  knowledge  of  some  of  the  statistical  properties  of  the  scattered 
waveform.  Specifically,  we  need  a  suitable  model  for  the  chaff 
cioud  from  which  we  can  deduce  various  statistical  properties  of  the 
received  waveform.  For  design  purposes,  the  most  important 
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information  may  be  considered  to  be  contained  in  the  received  signal 
correlation  functica  [233* 

* 

Consider  a  chaff  cloud  composed  of  a  large  number  of  dipoles 
which  are  moving  about  and  reflecting  energy  independently  of  one 
another.  In  addition,  each  dipole  is  assumed  to  be  continuously 
tumbling  (rotating)  about  an  axis  perpendicular  to  its  length.  Assume 
that  the  transmitted  signal  is  of  the  form 

(2. 37)  ' 

(2.38) 

is  the  complex  modulation  envelope  and  is  the  angular  carrier 

th 

frequency.  Let  a^(t)  he  the  reflection  coefficient  of  the  k  scatterer, 
and  r^  and  r^  be  the  distances  from  the  k^1  scatterer  to  the  trans* 

A 

mitter  and  receiver,  respectively.  The  complex  envelope  of  the 
received  signal  over  an  observation  interval  (-T,  T)  can  be  repre* 
sented  by  [24],  [25], 


fas  t 

eo(t)  *  Re[s(t)e  °  ] 

\ 

where 

if  At) 

•it)  =  |T(t)|e 


Zr.  its  ,.t 

+  ~>7(t-t  )e  dk  (2. 

V«-T-T> 


39) 


where  the  subscript  k  designates  the  k*^  scatterer,  t^  =  (r^  +  r^/c 


is  the  time  delay  between  the  transmitter  and  receiver,  *dk  Uth* 
Doppler  frequency  shift,  and  g^  is  a  complex  number  which  assy  be  a 

function  of  the  positional  coordinates,  the  antenna  beam  characterise 

*  '  * 

tics,  and  other  initial  parameters.  From  Eq.  (2.11),  the  reflection 
coefficient  can  be  written  as 

*kW  '  V  {5(?k-  V  *  ®<V  V  'x>*  t*  l2**rk**V3 

+  X. ( ?k.  \)  exp  [  -i2(*rkt  +ak5]|  (2.  40) 

The  following  simplifying  assumptions  are  made  regarding 
the  scatterers: 

(1)  The  scatterers  are  statistically  independent  and 
‘  identically  distributed. 

(2)  The  distribution  of  the  scatterers  obeys  Fo/.seoa 
statistics,  i.e.,  the.  joint  probability  for  n  reflections 
to  oceur  in  the  time  intervals  (tj,  tj  4  dtj), 

(t2.  t2  +  dt2) . (tn,  t#  +  dt^)  is  given  by 

P<V*2 . ln;n,dtId*2 . *dtn 

"  r  fT  , 

=  a  — j~  exp  [- J  v(t)  dt  ]  d^  (2.41) 

k*l  *■  -T  J 

where  v(t)  is  the  rate  of  signal  arrival  at  the  receiver. 


(3)  The  velocity  of  the  sca&terers  end  the.  magaltadt  of 
the  scattered  signal  are  weakly  correlated  within  the 
observation  interval  (*T,  T|u 

(4)  The  orientation  of  the  rotation  avis  is  not  greatly  affected 
by  the  motion  of  the  scatterers* 

( 5)  The  initial  orientation  angle  of  the  dipole  axis  la 
uniformly  distributed  in  the  interval  0  <ft^<  2tt, 


For  simplicity,  we  sssume  <z(t)>  *  0  and  the  correlation 


function  is  given  by 


s  2  <*M*1)*?t2)> 


-  T 


*Sc> 


.  *  -♦-»  -*-*-*-*  -t  •*  ^  W 

•J*  *  •  j  p(g.  »Q.  5,  *1.  a,  jud,  JUy)dg  •  •  *  djtr  lR*(tl)  s(t2)J 

(2.42) 


J•••J  p(g.ao,...,a3r)  dg... d»s 


+• 

=  j  ...Jp(i)  <**  ...J  ...Jp(*r)d»r 

•  dp  « 


with 
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(2.43) 
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viMr* 


K  *  \  <!•!*>  <1*  2> 


T  *  V*i 


(20 44) 
(2. 45) 


M(r)  =  C<c2>  +  <U2>*  (2T)+<L*>#  (-2t)] 

r  r 

to,  t 


(2.46) 


*  <t>  *  |  ?(»„) .  * 

Q  e  00 


=  characteristic  function  of  p(ts^)  (2.47) 


to  t 


*a  (»)  *  J  p(»r>  •  r  <*»t 


=  characteristic  function  of  p(i»  )  (2.  48) 

and  C,  U,  and  L  are  magnitudes  of  C,  U,  and  L,  respectively. 
Equation  (2.43)  expresses  the  lime-varying  correlation  of  the 
received  waveform.  However,  if  v(t)  varies  only  slowly  with  time 
and  if  the  observation  interval  is  long  compared  with  signal  duration, 
(2. 43)  reduces  to 


,+® 

R(t)  =  Kvo*b  (T)M(T)J  *;*(u-T)Mu)du  (2.49) 


where  is  the  average  rate  of  signal  arrival  at  the  receiver. 
Thus,  the  correlation  function  depends  only  on  the  time  difference 
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(tj-tj)  and  the  process  becomes  stationary.  Furthermore,  if  the 
dipoles  are  not  rotating  (»f  =  0),  (2.49)  reduces  to 

+** 

R(t)  =  Kv  *  (r)<h*>f  ^*(u-T)7(u)du  (Z,  SO) 

®  4  °  V* 

2  2  2  2 

where  <b^>  =  <C  +  U  +  L  >  is  the  relative  average  cross  section. 


2.  5  Evaluation  of  Scattered  Power  in  Spectral  Components 

The  computation  of  the  received  signal  correlation  function 
2  2  2 

requires  the  values  of  <  C  >,  <U  >,  and  <L  >.  These  quantities 
represent  the  relative  magnitudes  of  the  average  power  contained  in 
the  three  spectral  components.  Assume  that  all  orientations  of  die 
rotation  axis  of  the  dipole  are  equally  likely.  Then,  from  Equations 
(2.  !3)  through  (2. 19).  ws  find 

<c2>  =  ^  {i  +<l»/-  !?/)<iS;|2-  l^ftcos* 

*  4  co* 

(2.51) 

<U2>  «  <L*>  »  <S2> 

*  £  {'  Mil/- P/hPjI*-  1^1*)  «•» 

*  77  -  (VV'Vty  «•  Ml*}’ 


(2.52) 
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Th«  received  signal  depends  upon  the  receiving  polarisation 
as  well  as  the  transmitting  polarization.  Thus,  by  properly  choosing 
the  values  of  a^,  a^,  a^,  and  a^  as  in  Equations  (2.9)  and  (2.  10), 
the  scattered  power  for  any  polarization  combination  can  be  obtained. 
In  many  practical  situations,  either  linear  or  circular  polarization  ia 
employed.  For  the  case  of  linear  polarization,  we  have 

<C2>  *  ^  [1  ♦  7  (sin$o  sinqp^  +  cosqp^  cosqp^  cosp)2]  (2.53) 

<S2>  *  “r  [1  +  --(sin<p  sinqp*  +  cosqp  cos<p'  cos£)2]  (2.54) 

fU  J  ©  ©  o  o 

>  j  *2 

<bo>  *  [1  ♦  2  (sinqpQ  sinqp^  +  coscp^  cosqp^  eos$)  ]  (2.55) 

where  (p  and  qp'  are  defined  in  Fig.  1.  For  this  case,  it  is  noted  - 
o  o 

that  minimum  received  power  is  observed  when  ' 

(simp  simp*  +  cosqp  cosip'  cos£)  =  0  (2.56) 

o  o  o  o 

and  maximum  received  power  is  observed  when 

(simp  sincp'  +  cosqp  cosqp*  cos))  *  *  1  (2.57) 

O  O  O  O  “ 

Equation  (2.56)  can  be  satisfied  for  many  combinations  of  qp  ,  qp*  , 

o  o 

and  6,  while  Eq.  (2.  57)  can  be  satisfied  only  for  limited  values  of 

qp  ,  qp*  ,  and  6;  i.e. ,  (1)  6=0  and  the  transmitter  and  receiver 
o  o 

polarizations  are  parallel  to  each  other  and  (2)  MO  and  both  the 
transmitter  and  receiver  polarisations  are  perpendicular  to  the 
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plane  of  incidence  {©  =  ♦  90°  and  $>'  =  +  90°}.  Similar  conclusions 

o  ■“  o  — 

have  also  been  observed  by  Borison  [6]. 

For  the  case  of  circular  polarisation,  we  find 

<C2>  =  ~[<1  +  cos3)+  ^  (1  +  cosp)2]  (2.SS) 

<S2>  =  (l  +  cos&)2]  12.593 

<b2>  *  |  cos2*]  (2.60) 

where  the  upper  sign  corresponds  to  the  case  when  both  the  trans¬ 
mitter  and  receiver  polarizations  are  circularly  polarised  in  file 
same  sense  and  the  lower  sign  corresponds  to  the  case  of  the  oppe- 
site  sense.  The  received  power  is  a  function  of  the  bistatic  angle  3 
and  the  cross  section  is  independent  of  the  sense  of  polarisation. 

2.  6  Conclusions 

The  signal  scattered  from  a  rotating  dipole  consists  of  &re# 
spectral  components  -  a  carrier  component,  an  upper  sideband 
component,  and  a  lower  sideband  oompooent.  The  sideband 
components  are  separated  from  the  carrier  component  by  twice  the 
rotation  frequency  of  the  dipole.  The  amplitudes  of  the  sideband 
components  are  in  general  not  equal  and  hence  both  amplitude  and 
phase  modulation  of  the  scattered  signal  can  result. 


The  average  cross  section  is  not  affected  by  the  rotation  of 
the  diposes.  For  ail  polarisation  combinations  of  the  transmitter  and 
the  receiver,  the  average  cross  section  varies  between  the  limits  of 
77  0  and  7  a  ,  where  0  is  the  maximum  cross  section.  In  general. 

19  O  9  O  O 

the  polarization  of  the  transmitter,  the  receiver,  or  both  can  be 
varied  to  achieve  either  maximum  or  minimum  received  power  for 
a  given  bistatic  angle. 

The  correlation  function  of  the  received  signal  envelope  from 
a  collection  of  randomly  rotating  dij>ole  scatterers  depends  upon  both 
the  transmitting  and  receiving  polarizations,  the  transmitted  signal 
envelope,  and  the  statistical  properties  of  the  scatterers.  ■ 


Chapter  3 


A  MODEL  FOR  THE  RADAR  ECHO  FROM  A  RANDOM 
COLLECTION  OF  ROTATING  DIPOLE  SCATTERERS 

3*  1  Introduction 

In  radar  systems  it  is  often  necessary  to  detect  a  target  echo 
in  the  presence  of  other  unwanted  echoes  or  clutter.  Clutter  echoes 
include  signals  reflected  from  chaff,  surface  of  the  ground  and  sea. 
vegetation,  etc.  In  order  that  a  radar  receiver  can  be  designed  to  • 
operate  effectively  in  the  presence  of  clutter  interference,  it  is 
important  to  develop  a  suitable  theoretical  model  of  the  clutter. 
Siegert  [26],  Kerr  [13].  and  Lawson  and  Uhlenbeck  [27],  presented 
the  first  two  probability  distribution  functions  for  the  magnitude  of 
the  returned  echo  from  an  infinitely  dense  chaff  cloud.  Kelly  and 
Lerner  [24]  developed  a  theoretical  model  of  a  chaff  cloud  from 
which  they  derived  various  statistical  properties  of  the  returned  echo. 
The  signal  returned  from  the  chaff  cloud  was  considered  as  a  random 
process  and  the  scatterers  were  treated  as  points  with  variable  cross 
sections. 

The  basis  for  the  analysis  can  be  described  as  follows.  An 
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3* 


RF  pals*  is  transmitted  from  a  radar  toward  a  cloud  of  random 
scatterers  moving  about  and  reflecting  energy  independently  of  on* 
another.  In  addition,  the  cloud  is  assumed  to  have  an  overall  drift 
velocity.  Since  the  wavefront  intercepts  scatterers- at  different 
ranges,  the  echoes  returned  to  the  radar  will  arrive  at  a  rate  which 
depends  upon  the  local  density  of  the  cloud.  If  the  effects  of  multiple 
scattering  and  scatterer  rotation  are  neglected,  the  echo  signal  from 
a  particular  scatterer  can  be  regarded  as  a  Doppler -shifted  replica 
of  the  transmitted  waveform.  However,  any  change  in  the  orientation 
of  the  scatterers  can  cause  variations  in  the  returned  echo  power  and 
phase,  which  in  general  cannot  be  neglected.  It  is  the  purpose  of  the 
present  analysis  to  show  that  rotational  motion  of  the  scatterers  can 
have  significant  effects  upon  the  echo  waveform, 

"Chaff"  is  a  form  of  countermeasure  used  against  radar.  A 
chaff  cloud  usually  consists  of  a  large  number  of  metallic  strips 
(dipoles),  which  closely  resemble  an  assembly  of  random  scatterers. 
When  chaff  dipoles  are  dispensed  from  a  moving  craft  in  space,  the 
effects  of  ejection  forces,  body  instability,  and  other  aerodynamical 
properties  can  cause  the  dipoles  to  rotate.  Thus,  ignoring  the  effects 
of  scatterer  rotation  car.  result  in  an  incorrect  representation  of  the 
chaff  model,  although  in  many  instances  such  effects  can  be  tolerated 
without  greatly  affecting  the  results.  Ground  clutter  exhibits  echoing 


characteristics  similar  to  those  of  random  scatterers.  Rotational 
motion  of  the  scatterers,  or  equivalent  movements  of  branches, 
leaves,  grass,  etc.,  under  the  effects  of  wind  forces,  can  become  a 
major  contributing  factor  to  fluctuation  of  echo  intensity.  Karr  [13] 
and  Lawson  and  Uhlenbeck  [27]  presented  some  experimental  results 
on  the  measurements  of  chaff  models  and  ground  clutter.  In  the  case 
of  chaff  measurements,  it  was  reported  that,  in  addition  to  the  fast 
Doppler-beat  fluctuations  a  slow  secular  variation  was  almost  always 
observed.  The  origin  of  this  type  of  fluctuation  was  not  clear,  but  It 
was  conjectured  that  the  cause  could  be  due  to  the  rotation  of  dipoles. 
For  the  case  of  ground  clutter,  it  was  shown  that  the  width  of  the 
fluctuation  frequency  spectrum  was  not  proportional  to  radio 
frequency  as  would  be  expected  if  the  fluctuations  were  due  entirely 
to  the  velocity  distribution  of  the  scatterers.  At  some  wavelengths, 
the  fluctuation  spectra  approximately  resembled  those  of  randomly 
moving  scatterers,  while  at  other  wavelengths  the  measured  results 
differed  significantly  from  the  theoretically  expected  values,  Kerr 
indicated  that  this  anomalous  behavior  eculd  be  due  to  fluctuations  of 
two  sources  instead  of  one,  but  no  theoretical  explanation  was  offered 
to  clarify  this  phenomenon. 

A 

This  chapter  presents  a  refinement  of  the  theoretical  model 
for  the  radar  echo  from  a  random  collection  of  scatterers.  The 


scatterers  are  treated  ac  dipole*  instead  of  point*  with  variable 
cross  sections.  This  assumption  permits  the  determination  of  the 
effects  of  scatterer  rotation  which  have  previously  been  neglected 
£24]  -  [27].  The  theory  can  be  extended  to  the  study  of  other  classes 
of  clutter  signal  which  exhibit  similar  echoing  characteristics*  For 
example,  targets  making  up  ground  clutter  such  as  trees,  grass, 
etc. ,  will  move  and/or  oscillate  with  the  wind.  The  individual  blades 
or  stems  scatter  energy  like  lossy  dipoles.  For  our  present  work, 
we  ignore  the  loss  characteristics  of  the  scatterers.  The  relative 
magnitude  of  energy  reflected  by  individual  scatterers  is  taken  eare 
of  in  the  probability  distribution,  of  course. 

An  expression  for  the  correlation  function  is  derived  in  terms 
of  the  characteristics  of  the  transmitted  waveform,  polarization,  and 
the  distribution  of  the  scatterers.  Under  the  assumption  that  the 
cloud  density  is  slowly  varying,  the  process  becomes  stationary.  It 
is  shown  that  the  power  spectrum  of  the  returned  echo  consists  of 
three  closely  spaced  components,  one  centered  at  the  mean  Doppler 
frequency,  while  the  other  two  are  centered  at  the  sum  and  difference 
of  the  mean  Doppler  frequency  and  the  twice  the  mean  rotation 
frequency  of  the  dipoles.  The  fluctuating  characteristics  of  clutter 
echoes  are  also  determined,  and  the  theory  is  compared  with  some 
of  the  experimental  results  in  the  literature. 
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3.  2  Analysis  of  the  Model 

Let  the  coordinate  system  be  centered  at  the  radar.  The 
cloud  consists  of  a  random  collection  of  moving  dipoles  reflecting 
energy  independent  of  one  another.  In  addition,  each  dipole  is 
assumed  to  have  a  periodic  rotation  about  an  arbitrary  axis  perpen¬ 
dicular  to  its  length.  ^  The  polarisation  of  the  transmitted  wave  is 
assumed  to  be  uniform  throughout  the  illuminated  volume  and  the 
effects  of  multiple  scattering  are  neglected.  As  in  Equations  (2.  37) 
and  (2.  38),  assume  that  the  transmitted  signal  has  the  form,  * 

eo(t)  *  Re[s(t)e  °  ]  .  (3.1) 

as# 

where  s(t)  is  the  complex  modulation  envelope.  The  returned  echo 
from  a  collection  of  scatterers  located  at  (r^,  0^,  |^)(  k*l, Z  3,..., 
can  be  written  as 

x(t)  a  Re  [«(t)3 

Sc  ~ 

k 

.  exp  mdk)(t-tk)+ «(#k,  #k)))  (3.2 

*  The  term  "rotation"  used  here  implies  the  instantaneous  change  in 
orientation  of  the  dipole  with  respect  to  the  direction  of  the  radar 
during  the  time  interval  in  which  th-e  dipole  is  illuminated  by  the 
radar  beam.  For  the  purpose  of  calculation,  it  is  convenient  to 
represent  such  *  change  in  orientation  by  a  periodic  rotation  about 
an  appropriate  axis. 


49 


wktrt 


Kj  =  *  canitut 
C  *  power  gain  of  antenna^ 
s  reflection  coefficient 
2*o  ^ 

ir^  =  ^  ;fio  -  Doppler  frequency  s  lift 

=  overall  drift  radial  velocity  of  the  cloud 
v^  =  differential  radial  velocity  of  k  scatterer 
*  velocity  of  light 

=  phase  shift  depending  on  initial  parameters 

e  =  phase  shift  due  to  the  antenna  beam  phase 
characteristics. 

The  geometry  for  the  k**1  scatterer  is  sh-^wn  in  Figure  3. 
The  rotation  axis  of  the  dipole  is  described  by  the  angles  and  t^* 
where  5^  is  the  angle  between  the  rotation  axis  and  the  unit  vector 
e^  (*he  direction  of  propagation)  and  is  the  angle  between  the 
projection  of  the  rotation  axis  in  the  d  0-plane  (the  plane  perpen¬ 
dicular  to  the  direction  of  propagation)  and  the  unit  vector  e^>  The 
reflection  coefficient  takes  the  same  form  as  Eq.  (2.40),  i.e.  , 


If  the  antenna  has  different  gains  for  transmitting  and  receiving* 
i.e. ,  due  to  polarization  changes,  then  C  should  be  replaced  by 


*k(0  *  *ok  CC<?k.  r^)  +  UJ5k.  t^)  exp  C  +  i2(arkt+ak)) 
+  L<5k.  exp  [-i2(aykt  +  ftkU3 


(3.3) 


where  7  .  is  the  broadside  echo  voltage  of  the  lipole  with  the 
electric  field  vector  parallel  to  its  axis,  te^k  is  the  angular  rotation 
frequency  of  the  dipole,  ie  the  initial  angle  between  the  dipole 
axis  and  the  plane  perpendicular  to  the  direction  of  propagation,  and 

^  «v» 

the  quantities  C,  U,  and  L  are  functions  of  polarisation  and  the 
orientation  of  the  rotation  axis  of  the  dipole. 

If  we  assume  that  the  cloud  occupies  a  finite  volume  in  space 

and  the  mean  distance  between  the  cloud  and  the  radar  is  large  se 

that  the  difference  between  the  range  of  a  particular  scatterer  r^  and 

the  mean  range  r  is  small  compared  with  r  ,  we  can  ignore  the 
o  o 

positional  dependence  of  the  amplitude  of  the  echo  signal.  Hence, 

substituting  (3.  3)  into  (3.  2)  and  setting  r.  -  r  ,  we  obtain 

&  o 

M»  *  —  t  0(9  *  (7  7|t-t  !«pl+l<|>  (t-yj 

r  k«l 

o 

•  {c(§k»  +  U^k’  ^k*  *xp  ^  +  2 


+  L(§k*  \)  ®*P  C-iViCt-~)J  J 


(3.4) 


l 
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where 

«Pc(t)  *  (»0  +  *dk)  t  ♦  #k  ♦  «(0k,  #k)  (3.  5) 

•  » 

*  2(»rkt*«k>  (3. 4) 

The  same  statistical  assumptions  as  in  Section  2. 4, 

Chapter  2,  can  be  made  regarding  the  scatterers.  The  random 
values  aok#  6^,  0^,  ?k>  Pk»  &k»  and  i#^k  are  taken  from  tho 

joint  probability  distribution  where  it  is  presumed  that  a  scatterer 
is  present  at  range  rfe  and  time  *k/2.  It  is  shown  in  Appendix  n  that 
the  probability  distribution  of  the  echo  waveform  approaches  Gaussian 
as  the  rate  of  echo  arrivals  becomes  large. 


Following  the  same  procedure  as  in  Chapter  2.  the  received 
signal  correlation  function  is  found  to  be 


R(t)=  K  v  e  °  *  (t)M(t)  f  7*{u-T)7(u)du  (3.7) 

2  w  *d  J 


where 


K,“  -"4  |K  12<0*><|»  t*» 

2r 

o 


C3.#) 


v  =  average  echo  rate 
o 

M(T)  *  C<C2>  +  <U2>*  (2  t)  ♦<  L2>  •  {-2T)3 

w  '£.-; 

r  r 


CJ.9) 


il  ,» 


*  (’)  *  |  P(»4)»  *  «•, 


*  characteristic  function  of  p(e^) 

*•  (T|  *  I  *•%*•  *  4,t 


*  characteristic  function  of  p(a^) 


(3. 10) 


I3.1H 


and<C*’>,  <U‘>.  and<L*>  are  the  expected  values  of  |c|^,  |fj|\ 
and  |l|*.  respectively.  If  we  define  p{r)  *  R(t)/R(o)  where  R(e)  Is 
the  noise  power,  (3. 7)  becomes 


p(T)s  - --  A(t)*  (t)M(t) 

<b  *>  *d 

o 


(5. 12) 


where  <  b  ♦  U*  +  L*>  and 

P 


A(f) 


j  “s*(t  -t)‘s(»)  dt 

•  m 

J*  p«t)l#* 


(3.13) 


is  the  normalized  autocorrelation  ol  the  transmitted  signal  envelope. 
The  power  spectrum  of  the  echo  wayeform  can  be  obtained  by  taking 
the  Fourier  transform  of  (3, 12).  Thus, 


.  *  % 


4S 

g(u>)  *  — 7-  [<C2>F  (a-u»  )  +<U2>F  (u>  •  *  )  +<L2>FJ«-®  )) 
<b2>  e  ©  «  *  *  • 

°  (3.14) 

whtrt 

F.(J»)  *  J*  p(Jt^)  £(»-«d)d«d  (3.  IS) 

•  •  • 

+  « 

Fu(u>)  *  J  J  p{i>d)  p(*r)  £(tn-od- 2ar)  d*d<l«r  (J.  16) 

•  m  m  m 

Fj(aj)  *  J  J  p(*d)  p|*r)  f(w-*d+2(«r)  di»dd«r  (3.17) 

mm  mm 

*• 

and 

+• 

f  (uu)  a  J  A(T)*’**Tdt  (3.19) 


is  the  puwer  spectrum  of  t?vi  transmitted  signal  envelope.  The 

power  spectrum  is  symmetric  if  p(<*’d)  and  p(u»y)  are  symmetric. 

For  example,  if  p(x^)  and  p(<»^)  are  Gaussian  with  means  »d  and 

•**,  the  power  spectrum  can  he  interpreted  as  consisting  of  three 

closely  spaced  components,  one  centered  at  the  mean  Doppler 

frequency  (j;  +  ju  .)  while  the  other  two  are  centered  at  (is  +  is.  +  2<s  ) 
Co  o  r 

and  (u»  ♦  is .  -  2  'is  ). 
o  q  r 

It  is  noted  that  the  normalized  correlation  function  and  power 


spectrum  do  not  depend  on  polarisation  if  the  effects  of  scatterer 


rotation  are  neglected.  Thus,  if  is  3  0,  (3. 12)  and  (3. 14)  reduce  to 

9 

to  t 

PJT)  .  A(t)*  (T)e  °  (3.19) 

•  d 

*0(»)  *  (3. 10) 

These  results  are  equivalent  to  those  obtained  by  others  [24],  [29], 
provided  that  the  echo  rate  is  nearly  uniform. 

The  average  backscattered  power  contained  in  the  three 
spectral  components  can  b?  e  valuated  by  using  Equations  (2.  51) 
through  (2.60).  Since  the  scattered  power  is  polarisation  dependent, 
the  results  for  the  following  special  cases  are  summarized  below: 

(1)  linear  transmit  •  linear  receive 

(2)  circular  transmit  -  circular  receive 

(3)  orthogonal*  linear  transmit  and  receive* 

(4)  orthogonal  circular  transmit  and  receive 


*  This  implies  that  the  transmitter  and-  receiver  polarization  are 
orthogonal  to  each  other,  i.  e. ,  e^  •  •*  *  0. 
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Case 

<c*> 

<u2> 

<l2> 

A 

•r» 

V 

(1) 

2 

1 

1 

1 

13 

30 

30 

5 

(2) 

1 

30 

1 

20 

1 

20 

2 

IS 

(3) 

1 

1 

1 

1 

60 

40 

40 

IS 

(4) 

7 

1 

1 

2 

60 

120 

120 

IS 

CS.21I 


Other  cases  can  be  evaluated  in  a  similar  fashion, 

3.  3  Characteristics  of  Echo  Fluctuations 

The  echo  power  fluctuates  with  time  because  of  the  motion  of 
the  scatterers.  For  simplicity,  we  assume  in  the  following  analysis 
that  the  mean  number  of  scatterers  at  a  given  range  is  large  so  that 
the  process  becomes  essentially  Gaussian.  The  first  and  second 
probability  distribution  for  the  echo  power  are  thus  given  by  [13], 
[26]  -  [27X, 


o  d  P 

Wj(P)  dP  *  s  ~ 

o 


(3. 22) 


P1  +  P2 


P0d-P'2)  2p*/CT?  dPl4P2 

w,<pi.p,.T)dP  dp  *  e  w  i0( - T-J~z - T 

2,2  12  °VPA(1-P2)  'pAi-P'*) 


(3. 21) 


"rsMj&a  ' 

'  *-  -3£r;<' 

\  r  -T-.  ./•  .■ 
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*?here  IQ  is  the  seroth  order  modified  Bessel  function  of  the  first 
kind,  P0  is  the  average  echo  power,  P^  =  P(tt),  Pj  *  P(t^|, 

T  *  ^*tj.  »»d  p*  is  the  correlation  coefficient  between  the  echoes 
of  successive  pulses.  Prom  (3. 12),  we  obtain 


P'(t) 


1 


<»,*> 


ifll  T 

Re[s  °*  (t)M(T)) 


13.21) 


The  correlation  function  for  P  can  be  defined  as 


I(T) 


<  p.  *  P«HP*-  P«)>  <P.  P,>  -  P, 


1  0*2  0 
*  PI2>  -  p« 


1  2  0 

<pi2>-p« 


(3.25) 


Using  (3. 20)  and  (3.  21),  we  obtain 


<lV*  Jwi'pi*p,l,pi’  2P« 

0 


<pipi> *  If  vy'w>*,,pz<*vp2*  po«,+",2» 

0  0 


Thus,  (3. 25)  becomes 


I(t)  *  Cp’lt)]' 


(3.26) 


Now,  we  assume  that  both  the  Doppler  shift  and  the  rotation  frequency 

have  a  Gaussian  probability  distribution  with  means  i.  and  •  B 

®  If 
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respectively.  Thus. 


p(>»,)  - J exp  r.  («  -»  )*/2o  *  1  (3.  27) 

d  /TiTo  1  d  d  d  J 


p^j  *  — ■ —  «*p  r-(* -»>w  i 

r  /TS*«  L  *“  r  r  J 


/Tit*  o 


(3. 28) 


From  (3. 10)  st>d  (3.  11).  the  characteristic  functions  of  p(is  )  and 

d 


pja^j  become* 


1_  <T> 
•a 


—  ,  z  z 

‘V  -*°aT 

e  e 


13.29) 


to  t  -|o2t2 
(t)  -  •  r  .  r 
r 


(3.30) 


Substituting  (3.9).  (3.  29)  and  (3.  30)  in  (3.  24)  and  suppressing  the 

‘<V"d,T 

phase  factor  e  (which  corresponds  to  mixing  the  returned 

A 

signal  with  an  oscillator  of  frequency  JU^+ajj),  then  (3.  26)  becomes 


I(T)  *  « 


T  <C*>  .  '2°r  '  f< .  <l?» 

<b2>  .^<b*>  <b*> 

o  o  o 


2  %  2 

Since  <t)  >  *  <L  >  =  < S  >,  (3.  31)  can  be  written  ns 


(3,  31) 
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I(t)*o 


2  2 

-v 


2  2  -2o2t* 

<C  >  ,  ,  <S‘>  r 
- r —  +  2  — ~  e  cos  2«s  t 

<b2>  <b2>  ' 

o  ® 


(3.  32) 


The  frequency  spectrum  of  the  power  fluctuations  can  be  obtained  by 
taking  the  Fourier  transform  of  1(t)„  Tims, 
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where 
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c2 *  o?  +  4a 2 
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*.  «.f  *  *•“ 

b  d  r 


(3. 33) 

(3.34) 
(3.  35) 
(3.36) 


One  of  our  objectives  was  to  clarify  theoretically  some  of  the 
anomalous  phenomena  observed  during  early  experimental  studies  of 


SI 


clutter  [13].  Figure  4  illustrate*  the  correlation  function  lor 
wooded  ground  clutter  measured  at  9.  2  cm  and  at  wind  specie  of 
SO  mph.  The  experimental  data  were  taken  from  Kerr  [IS.  p,  387], 
The  theoretical  values  are  shown  by  the  dashed  curve,  which  can  bn 

obtained  from  (3.  32)  and  (3.  21)  with  <C2>/<b2>  *  |  and 

o  3 

2,21 

<S  >/<b  >  =  7  for  linear  polarisation.  Thus,  we  have 
o  o 


-0.0289  x  JO4?1 


-0.  309  *  I04?2 

e  cos 


(3.37) 


This  corresponds  to  o.  =  17  rad/s  and  m  *  0  *39.3  rad/s.  The 

o  r  r 

corresponding  frequency  spectrum  of  the  fluctuations  can  be  obtained 
by  taking  the  Fourier  transform  of  the  correlation  function,  and  is 
shown  graphically  in  Figure  5. 


It  should  be  mentioned  that  the  measured  correlation  functions 
were  obtained  by  radar  envelope  measurements.  However,  Kerr 
[13]  has  indicated  that  the  difference  between  the  correlation  func¬ 
tions  for  envelope  and  power  is  small,  the  maximum  deviation  being 
0.  027.  Hence,  by  neglecting  this  small  difference,  the  two  correla¬ 
tions  can  be  considered  equivalent. 


It  is  noted  that  the  curves  deviate  significantly  from  the  n*ugl 
expected  Gaussian  shape.  Instead,  they  drop  very  slowly  toward 
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9.  2  cm  (Heavily  Wooded  Terrain  at  Wind  Speeds 
of  50  mph) 
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aero.  The  anomalous  appearance  of  the  correlation  function  is  due 
predominantly  to  the  oscillating  character  of  the  scatterer*.  K«r?p 
has  pointed  out  that  a  small  positive  asymptote  has  been  subtracted 
from  the  original  measured  correlation  function  in  an  attempt  to 
correct  for  the  slow  secular  variations.  The  correlation  function 
has  been  renormalized  to  a  value  of  1. 0  at  T  =  0.  Hence,  an  exact 
comparison  between  the  calculated  and  measured  values  cannot  be 
made.  However,  the  close  resemblance  between  the  shapes  of  the 
theoretical  and  experimental  curves  indicates  good  agreement  of 
theory  with  experiment. 

3.4  Conclusions  , 

When  a  radar  operates  in  an  environment  of  clutter,  it  is 
important  to  determine  all  possible  effects  which  can  limit  the 
effectiveness  of  the  receiver.  If  the  echo  waveform  were  precisely 
known,  it  would  be  possible,  in  principle,  to  synthesize  a  filter  which 
optimizes  the  signal-to -noise  ratio  at  the  output  of  the  receiver.  In 

this  chapter,  we  have  emphasized  the  effects  <?f  scatterer  rotation. 

t 

i  * 

The  following  conclusions  can  be  drawn. 

1)  If  the  spectral  spread  due  to  the  Doppler  velocity 

fluctuations  of  the  scattered  is  large,  then  the  effects  due 
to  scatterer  rotation  can  he  neglected.  Receiver 


ss 

characteristics  will  be  primarily  determined  by  the 
Doppler  spectrum  of  the  clutter* 

2)  If  the  Doppler  spread  is  small,  then  the  contributions  due 
to  rotations  of  the  scatterers  can  become  significant. 

This  was  demonstrated  in  the  example  considered  at  the 
end  of  Section  3.  3. 


.....  -  -  - - 
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ON  SCATTER  COMM  UNICATION  VIA  RANDOM 
TIME-  VARIANT  "DISCRETE  SCATTER ER"  CHANNEL 


4.  1  Introduction 


In  m.'ny  radar  applications,  communication  may  be  estate 
lished  between  two  distant  points  by  means  of  electromagnetic 
scattering  from  a  collection  of  scatterers  randomly  distributed  la  a 
region  of  space.  An  elementary  scatter  communication  system  Is 
shown  in  Fig.  6.  The  region  of  space  occupied  by  the  collection  of 
scatterers  is  referred  to  as  the  scattering  medium.  An  example  of 
such  a  medium  is  a  cloud  of  resonant  dipoles  (chaff)  C7]-[ll3«  Tbs 
basic  equation  relating  the  system  parameters  is  the  bistatlc  radar 
equation. 


cT  gr  r 

~3  FT 

<4">  rTr* 


where  the  subscripts  R  and  T  refer  to  the  receiver  and  transmitter 
parameters,  respectively,  and  the  quantity  a  is  the  scattering  func¬ 
tion  which  depends  upon  the  characteristics  of  the  mediron. 
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Figure  6. 
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Scatter  Communication  System 
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b  general  the  scattering  mediun  variee  with  both  posit*  ca 
and  time.  When  viewed  from  the  transmitting  and  receiving  terminals, 
the  scatter  channel  acts  like  a  linear;  random,  time-varying  four* 
terminal  network.  Our  objective  is  to  develop  a  suitable  response 
function  for  this  network. 

In  the  analyses  of  random  scattering,  most  investigators 

a 

Ignore  the  time  dependence  and  are  concerned  only  with  the  mean 
power  which  may  be  provided  by  the  channel  with  a  sine-wave  excita¬ 
tion,  For  example,  in  Eq.  (4.1),  a  usually  represents  the  averags 
cross  section  of  the  scattering  medium.  In  a  practical  situation,  the 

medium  it  time- varying,  the  transmitted  waveform  is  a  modulated 

« 

carrier,  and  the  received  signal  is  passed  through  various  filters. 

The  design  of  such  filters  requires  more  information  than  the  mean 
channel  power  outgo*  or  ih®  average  cross  section  of  the  scattering 
m«4Uw-  . 

The  present  analysis  is  an  extension  to  that  of  Kelly  [  141. 

% 

The  general  theory  of  electromagnetic  scattering  from  a  random 
collection  of  small  metallic  scatter**®  is  developed.  The  antenna 
gals  functions  are  included  in  the  derivation  of  the  channel- function 
in  order  to  take  explicit  account  of  the  possibility  that  scatterers  may 
flew  in  and  out  of  the  volume  illuminated  by  *he  two  antenna  beams. 
Arbitrary  transmitting  and  receiving  polarisation.*  arc  considered. 


and  the  effects  of  scatterer  rotation  are  also  taken  Into  account* 
Finally*  specific  results  are  obtained  for  the  case  of  dipole  scattering* 


4.2  Derivation  of  the  Scattered  Field 


Consider  the  scattering  from  a  collection  of  metallic  scatter* 
ere  which  are  contained  within  a  finite  volume  V.  Assume  that  the 
origin  of  coordinates  is  also  contained  within  V  and  that  all  sources 
which  produce  the  primary  field  are  external  to  this  volume.  The 
fields  everywhere  satisfy  the  Maxwell's  equations* 


V  X  E  *  -4 

o  ot 

—  dE  — 

, ,  H  .  ,o  w  ♦  j 


(«.*» 


where  J  is  the  source  current  which  gives  rise  to  the  primary  field* 
If  the  scatterers  are  perfectly  conducting*  the  boundary  conditions 
n*H  =  0  and  n  x  E  =  0  are  satisfied  on  the  surface  of  each  scatterer* 


It  follows  from  (4.2)  that  E  satisfies  the  vector  Helmholts  equation* 


-  a2E 

V  x  V  x  E  +  u  c  *  -w 

o  o  at2  ® 


(4.3> 


We  assume  that  the  tir^e  dependence  of  the  fields  and  Its 


sources  can  be  resolved  into  Fourier  components* 


J 
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Equation  (4.3)  becomes 


V  X  V  X  E  -  k2E  *  -l®  ji  J 

o 

where  E  =  E(r,u>),  J  =  J(r,u>),  and  k^  =  ®^ii  «  • 

o  o 


(4.5) 


Application  of  the  vector  Green's  theorem,  together  with  the 
proper  boundary  conditions  and  the  radiation  condition  for  r  **«•,  to 
(4.  5)  yields  (Appendix  III) 


E(r)  *  "(7,7')  •  J(?)dV 


where 


+  f  r(7.7")  •  n  X  H(7")  dS»  } 

JS  J 

-tv|7-7l 


(4.6) 


"<r*%>  =  cr+ivv)e— . - 


4tv  J  r-r  I 
'  o 


is  the  dyadic  free  space  Green's  function  [28]  and  I  is  the  identity 

dyadic.  Other  field  quantities  are  related  by  Maxwell's  equations. 

The  first  term  on  the  right  side  of  (4.6)  is  the  primary  field,  E  (r). 

o 


Hence,  the  scattered  field  is  given  by 

_  _  r  — _  _ 

Es(r)  =  I7c~J  r<r'r">'  5  *H(r")dS" 
o  S 


(4.8) 


where  the  integration  is  carried  over  the  surface  of  the  scatterer. 


The  geometry  for  a  single  scattcrer  is  shown  in  Fig.  7.  W* 


Figure  7. 


Coordinate  System  for  a  Single  Scatterer 


■■  **i#t*1 
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8*4  x*  fe3  Jhe  v*c lor  distance  from  the  coordinate  origin  to  the  center 
of  the  scatter*?,  and  r"  be  the  vector  distance  from  the  coordinate 
origin  so  a  point  on  the  surface  of  the  scatterer.  In  addition,  we  let 
r  »  r**«.r's  Then,  in  the  far  *one,  j  r-r"|*«  |7-7‘ !  *.£'•?•  *h«*re 
&*  is  a  unit  vector  in  the  (r-rv)  4irectic«it  rnd 


*.  -  ~  -lk|.?«r'|  s  +ik£' 

r{r,r'')  M  — - (I  -  a'a')  *  * 

4r?7«7*J 


•  r 


(4.9) 


Equation  (4.8)  becomes 


E.<?)  - 


-Ik  |  r-r' 1 
e  1 

4tt  17-7'  I 


A,  A. 

a'xa'  x 


{  [  n  x  H(7")1e+lh  *9’  r»  dS"  (4. 10) 
JS 


If  the  effects  of  multiple  scattering  due  to  the  presence  of  other 
scatterers  are  neglected  and  if  the  3catterer  dimension  is  small  com¬ 
pared  to  the  wavelength,  the  surface  current  induced  on  the  scatterer 
may  be  approximated  by  that  due  to  the  primary  field  aloe*.  If 
k|  r#  |  is  sufficiently  small,  the  exponential  inside  the  integral  may. 
be  replaced  by  il  +  ik  a'  •  r^).  With  these  approximations,  equation 
(4. 10)  becomes 


Vr> 


.2  -Ik!  r-r' I 

k  e _  __ 

*o  •  4n  lr-r*| 


a’xa'n  f  [  n  x  H  (7")]  (1  ♦  ik  a' .  7  )  dfS” 
J.  o  s 


(4.11) 


Next  it  can  be  shown  that,  [29} 


f  (a'.r  )[nx  H  {r")]dS"  *  -  5-a*  x  f7  x[nxil  (r”)]dSH 

Jg  •  O  Z  ,L  S  © 


(4. 12) 


Hence 


.2  -Iklr-p*  I 
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r  _  W  - 

- -  a'x(a'xp)  +  -a'Xm  1  (4.  IS) 

I  L  • 


«  I  r-r*  l 

where  p  and  m  are  the  induced  eleetric  and  magnetic  dipole  moments, 
respectively: 


P  *  J  CnxHo(7»)}dSM 

S 

m  =  i  J  7g  x  [nxHo(7")}dS" 


(4.  14) 

<4.151 


Thus,  in  the  far  zone,  the  scatterer  acts  as  though  It  were  both  an 
electric  and  a  magnetic  dipole  in  a  uniform  static  field.  In  classical 
electromagnetism  the  induced  electric  and  magnetic  dipole  moments 
are  also  expressed  as  linear  vector  functions  of  the  applied  fields 
r1)  and  H#(r'),  respectively  [30]-  [31],  as  follows: 


P  s  «ea*  Eo(y,) 


m*  *  *  H#<r‘) 
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where  a  end  8  are  the  dyadic  electric  and  magnetic  polarizabilities 
which  depend  on  the  size,  shape,  and  orientation  of  the  scattetfer* 


Now  assume  that  the  scatterers  are  located  at  points  rl 

1 

where  i  *  1,2,..., n.  By  Eq.  (4. 131  with  p  and  m  given  by  Equations 
(4. 16)  and  {4. 17),  respectively,  the  total  scattered  field  becomes 


.2  n  -iklr-r'f 

E.<7>  -  *  ! - L 

4  i*i  I  r-7: 1 


{  fj  X  Cij  x  ^•Eo(7j)l  +  .JxB^  Ho(rj|)  H.  l*» 


In  the  far  zone,  r  is  far  from  V  so  we  may  put  a'j  *  a^  and  I  7-7j|  * 
r  -  r  j.  With  this  substitution,  (4. 18)  becomes 


__  2  «  s  c  r  l 

E  (r,«)  *  --2-;  - -  Z  « 

51  4«c  r  i*l 


{  %xCarxSl.Eo<7j,»)]+~.arxIi.I7o(7j,»)}  (4.19) 

o 


2  2  2 

where  we  have  replaced  k  by  e  /c  . 


'  In  order  to  put  (4. 19)  in  th?  form  of  a  stochastic  integral, 
use  the  method  of  E.  J.  Kelly  >  14l.  Kelly  assumes  that  the  position, 
size,  shape,  and  orientation  of  any  scatterer  can  be  represented  by 
a  set  of  parameters,  x^x^, . . .  ,x^,  which  we  take  to  be  the  eoordi- 
nates  of  a  point  x  in  a  v  -dimensional  space  X.  This  space  is  a 
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direct  product  of  the  ordinary  three-space  V,  which  describes  tbs 
position  of  the  scatterers,  and  a  parameter  space  Y,  whose  coordi¬ 
nates  fix  the  size,  shape,  and  orientation  of  the  scatterer*.  Tkea, 
specification  of  a  point  x  in  the  one  scatterer  phase  space  is  reprs- 

^  -»  «S  •# 

sented  by  x  =  (r,y),  where  r  locates  the  scatterer  and  y  describes 

-*  III 

all  other  parameters*  If  locates  the  i  scatterer  in  its  phase  • 
space,  one  writes  =  afx^,  ^  =  8(xj),  and  7*  *  where 

—  -•  *•  *.  a 

a(x )  and  0(x )  depend  on  all  but  the  spatial  components  of  x  and  r'(x) 

has  components  equal  to  the  positional  coordinates  included  Inx. 

Next  we  introduce  an  integral -valued  measure,  N(S),  of  measurable 
sets  in  the  v  -dimensional  phase  space  X,  which  describes  the  collec¬ 
tion  of  scatterer*.  If  S  is  a  measurable  set  in  X,  we  take  NfS)  to  be 
the  number  of  scatterers  for  which  belongs  to  S.  In  this  .wanner, 
Eq.  (4. 19)  may  be  written  in  the  form  of  the  Stieltje's  integral. 


£(  r ,  uj)  *  -  - g 

4nc 


-l  -  r 
Z  c 

r*  e 


—  l 


+  i—  £  •  r!|x) 
c  r  1  ' 


£  Xft(7> .  El  r'(x))]  +  ax  9(7)  •  H  (7'(7). •)}  dN  £) 
l  r  r  o  «^c  r  o  l 

(4.  20) 

Eq.  (4.4)  or  taking  the  Fourier  transform,  we  find 
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4.  S  Derivation  of  the  Channel  Function 

Attume  that  a  transmitter  is  located  at  a  vector  distance  r? 
from  the  origin  and  that  the  collection  of  scattered.  which  produce, 
the  scattering,  is  in  the  far  sone  of  the  transmitting  antenna.  l£  a 
real  signal  f(t)  is  applied  to  the  input  terminals  of  the  transmitter, 
then,  in  the  scattering  region,  the  transmitted  Held,  can  be  repre¬ 
sented  by 

•  * 

Ejr.  t)  -  kJctM  •V["lVt",T'] 

H0(r. «)  •  f  V['4V  ‘'’Vi 

0 


>■ 


I 

I 

X 

i 

$ 


where 


ti  .  In  j  c  is  the  intrinsic  impedance  of  free  apace*  la  the  gal. 
o  7  o  o  _ 

function  of  the  tr.namltting  antenna.  *T  i.  a  unit  .ector  in  the  e T 
direction,  l?  ia  the  electric  polarisation  .ector,  and  ^  *  -i,*  *T 
(the  «».«  1.  propagating  in  the  -i,  direction).  Suhatitutln,  (t.t*)  » 


(4.22).  we  obtain 


4 


*9 


E§(r*0  *  j  I  CGjC'  U))]^2  W(x)£  [  t-s- 
X  — ^ 


6f,*Vc  (W‘*,{*nd*  dN<**k> 


C4.29J 


where  6{>)  it  the  Dirac  delta  function.  In  this  form  represent*  a 
superposition  of  plane  waves,  propagating  in  the  direction  of  the 
receiving  antenna,  with  slowly  varying  amplitudes.  Let  e  be  the 
polarization  vector  of  the  receiving  antenna  and  be  its  gain  foac>< 
tion.  Then,  by  flux  considerations,  the  plane  wave 


propagating  over  the  receiving  antenna,  will  produce  a  signal 


/5=Vgtgr  «V  "'«<*> 


at  the  antenna  terminals,  where  G  is  the  gain  of  the  receiving  awten- 

K 

na  in  the  -a^  direction.  Since  the  scatterers  are  in  motion,  the 
scattered  signal  will  in  general  have  a  time-varying  amplitude  and 
phase.  If  we  assume  that  the  receiving  antenna  is  sufficiently  wide¬ 
band  to  follow  these  variations,  the  total  received  signal  will  be 


h(t)  *  ^===“  J^j  CGT(7'(x))GR(7*(x))31/2CeK.  w«r>l  f*( t-e) 


ftC*-»o+ c  ^T+  da  «***(«•*& 


(4.  JO) 


For  convenience,  we  let 


-  1/2  i/2  - 

[CT(r')CR(r'n  »  C0otCor5  G(r-| 


(4.31) 


where  end  ere  the  maximum  gain  of  the  transmitting  and 
receiving  antennas,  respectively,  and  G(r')  is  the  normalised  gala 
function  relating  the  product  of  the  two  antenna  gains.  Substituting 
(4.28)  end  (4.31)  hi  (4.30)  and  using  the  value  of  K#  from  (4.25),  we 


h(t)  =  —=  i  ?  C(7'W\[eJi"Hx)- P(*)‘  Syl 
V  4nc2  "X  -• 


f  (t*~s)  ^  <*T+ ®R)*  *'(*)]  ds  dN(x,t) 


(4.32) 


where 


gotgorx 


(4.33) 


and  we  have  made  use  of  the  fact  that  eR*  aR=  0  and  hR*  eR. 

If  the  overall  scattering  process  is  regarded  as  a  random 
time-varying  channel,  we  may  put 
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This  is  a  generalisation  of  a  similar  result  obtained  by  Xelly  [14]. 

In  (4.37)  the  collection  of  scatterers  behaves  as  it  were  a  dispersive 

medium.  The  response  of  the  medium  is  described  by  the  channel 

function.  In  the  stationary  case  N  is  independent  of  time,  and  the 

2 

frequency  response  is  proportional  to  s  since  our  analysis  has  been 


based  on  the  static  polarizabilities  of  the  scattered.  If  the  scatter* 
ers  are  not  small,  our  procedure  remains  valid  if  «e  allow  the 
polarizabilities  to  depend  on  frequency.  In  this  case,  Eq.  (4.37) 

m  •+  **  -*  m  m  •* 

still  holds  with  a(x)  and  B(x)  replaced  by  S(s,x)  and  B(s,x). 


4. 4  Example:  Dioole  Scattering 


Consider  the  scattering  from  a  collection  of  axially  symmetric 
scatt erers  (i.  e.  ,  cylindrical  dipoles).  The  electrie  and  magnetic 
polarizability  dyadics  for  axially  symmetric  scatterers  can  be  reprs* 
seated  by 


a  *  d  d  ♦  p  p 


B  *  ♦  B2pI 


(4.39) 


where  d  and  P  are  unit  vectors  parallel  and  perpendicular  to  the- axis 
of  the  scatterer,  respectively.  The  components  of  the  dyadics  a  and 
B  are  complex  function  of  frequency  and  the  size  and  shape  of  ths 
scatterer.  The  transmitting  and  receiving  polarization  vector#  can 
be  written  as 


«T  s  (3*  ?T)i  +  (?•  $T)p 
*R  =  (3*  eR)d  +  (P*  SR)p 


(4.40V 


Hence 


V  *  *  *x=  t<ai-a2H5*  V  +  V  V  *Tn  (4-4l) 

and 

V  B- V  nBj-fl^Hd.^Hd.hT*  +  tyV  VI  (4.42) 


For  thin  dipoles,  the  magnitudes  of  <1.,  0  ,  and  0  are  small 

2  A  2 

compared  with  a^.  Then  . 

-  •*«. 

a<«,*)«  =-7  «T>  <4,«» 

V4fl  C 


and  the  channel  function  becomes 

V* 


K(t,«)  = 


V4n  c 


3 


J. 


C °,|(®t*Hd  •  eR|(5-  eT)]  dN(x,t) 


(4.44) 


where  a^(w,x)  depends  on  all  but  the  orientation  'oordinates  included 
«• 

In  x. 


Rotational  motion  of  the  scanners  can  be  taken  into  account 

* 

by  allowing  the  orientation  vector  d  to  become  time-dependent.  It  has 

been  shown  .n  Chapter  2  that  the  quantity  {3.  e  JlJ  -  S  )  may  be  put 

ft  T 

•  » 

in  the  form 


(i.  sTi  *  is*  sr l.-“*  ) 


<4.«» 
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where  C,  U,  and  L  are  functions  of  the  orientation  of  the  rotation 
axis  of  the  dipole,  and  f  is  the  instantaneous  orientation  angle  of 
the  dipole  axis.  Comparison  between  (2.28)  and  (4.43)  shows  that 
the  quantity  j  Q(o„x)|  has  the  significance  of  the  cross  section  of  an 
individtr&l  dipole.  If  all  orientations  of  the  rotation  axis  are  equally 
likely,  the  average  cross  section  is  given  by 

<  |  Q{»,x)  j2  >  *  0  <  l.2  >  (4.46) 

8  O 

where  0  is  Che  maximum  cross  section,  and 

o 

<  b2  >  =  <C2+  U2*  L2  > 

O' 

which  has  been  evaluated  in  Chapter  2  (Eq.  2.35).  In  the  preseut 
chapter,  the  bistatic  angle  6  (the  angle  between  the  directions  of  tho 
transmitter  and  receiver,  viewed  from  the  scattering  region)  is 

given  by 

cos  0  *  #  •  *  (4.47) 

T  1C 

4.  5  Statistical  Properties 

In  the  previous  sections,,  we  have  derived  results  for  ths 
scattering  f  rom  a  definite  set  of  scattercrs  in  space.  Ir.  this  section, 
we  will  examine  some  of  the  statistical  properties  of  the  channel 


output. 
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The  collection  o l  scattered*  ii  described  by  an  integral* valued 
measure,  N(S,t),  of  seUi  in  the  one-scatterer  phase  space  X.  Since 
the  scattering  is  random  in  nature,  N(5,t)  becomes  a  randomly  time* 
varying  random  measure  of  sets  in  X.  Poisson  statistics  are 
assumed  for  N(S,t)  so  that  the  random  variables  N(S^,  t),. . .  ,  N(S^»t) 
are  independent  if  the  sets  , . . . , are  disjoint,  and  that  the 
probability  that  N(S,t)  equals  the  integer  n  Is  given  by 


Pfr) 


[  <N(S,  t)  >T“ 


exp  [  *■  <  N(S,t)  >  3 


(4.40) 


where  <  N{S,  t)  >  =  E  N(S,t)  denotes  the  average  number  of  scatter* 
ers  in  set  S  at  time  t.  We  assume  that  the  measure  <  N(S,t)>  has  a 

N  • 

S« 

density  <n(c,t)>  so  that 


<  N(S,t)  >  *  j  <  n(x,t)>  dX 


(4.491 


where  dX  ia  an  infinitesimal  element  of  volume  in  the  phase  space  X. 
If  we  let 


■  then 


Wt  *  j  Wj(x)  dN(x,  t)  ,  1*1,2. 


E  Wf  =  j  Wj(#j  <  n(x,  t)  >  dX 


(4.50J 


U  51) 


g 

The  symbol  E  also  stands  for  the  expected  value  or  ensemble 
average.  We  will  use  these  two  notations  interchangeably. 


V.rwj.j  W*(3<n(?.D> 


flX 
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{4.  52) 


Cov  Cwj.^1  *  J  Wa(x)W  (x)  <n«x.t)>  dX 


(4,53) 


Now,  from  (4. 34)  *nd  (4. 36),  the  channel  response  June  Pice 
may  be  expressed  in  tw o  useful  forms, 

Mt)  *  J  k(t,  s)  f (t-s)  ds 


1  t* 


let 


2^  j  K(t,e)F(e)  e  do 


where 


k(t#,)  *  h  l  *<*••>  •  "*  «*« 


we  find 


2  MO  =  j  F  k{t,  s)  f(t-s)  ds 
o 


1  f  Ik* 

2n  J  E  K(t,o)F{e)  e*""  de 


and 


<4.  54) 


»4r  55) 


(4,  56) 


Cov  [  h(t),h{t')1  *  J  j  Cov  [  k(t, .),  k(f.  .)1f(|..,f(t..,.)  d.  da* 

»  o 

j  ♦" 

'  Z  J  J  C°v[K(t,u»)#K*(t,,«»)]ei^6"*  1  V(o)F*(e,>  dode* 
••  •• 


(4.57) 


The  two  channel  covariances  are  i elated  by 
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Cov  [k(t,  a),  Mt', •')]  * - 1  J  jc©vtK(t,n),K*Ct',»'))el<"*”*  *  \l«d©* 

3fl)  .«  .*» 

(*.») 


Equation  (4.  56)  implies  that  the  mean  output  to  f(t)  Is  d*ttf* 
mined  by  the  mear  properties  of  the  channel.  From  (4.37)  and  (4*  51), 
we  obtain 


E  K(t,  •)  * 


A  • 

.  o 


-lot. 


I 


i  = 
c 


<W 


r'(x) 


G(r'(x))  Q(<u,xJ  <  n(x,  t)  >  dX  (4.59) 


Recall  that  x  =  (r,y)  so  that  dX  -  dV  dY,  and  Q(»,x)  depends  only  on 
the  y-parameters.  Furthermore,  in  many  situations,  it  is  possiblo 
to  put 

<  n( x,  t)  >  =  f(y)  <  n(r,t)  »  (4.60) 


where  <  n(r,t)  >  is  the  r'ean  number  density  in  the  ordinary  space? 
and  f(y)  is  the  probability  density  on  the  y-parameters,  which  may 
also  be  time-dependent.  Then,  Equation  (4.  59)  can  be  written  in  the 
form 


E  K(t,«)  =  Afte  °  <  Q(ei,y)  > 


i-  (a  +  a  )»  7‘ 
(•  e  T  R 


•v 


G( 7‘)  <  n(7',t) >  dV*  (4. 61) 


where 


In  the  dipole  scattering  example  of  Section  4.4,  If  we  asetsma 


a  uniform  probability  of  all  orientations  at  any  one  time,  we  find 


<(d.eR)(d.ST»  *  j(«R*STi 


<  Q(®»y)  >  =  j  O 

3  Vine2  1  R  T 


(4.M) 


where  eR  and  e^  must  refer  to  the  same  coordinate  system.  In  term* 
of  the  notations  used  in  Section  2.2  of  Chapter  2,  we  put 


e  -*  a  A  ♦  a .  ■£ 
T  r  I 


e  _  i  1 

e_  =  a '  *  a',  •** 

R  r  I 


Here  e,^  and  eR  are  the  polarization  vectors  referred  to  the  trans* 
mittsr  and  receiver  coordinate  systems,  respectively.  With  the  aid 
of  the  rotational  transformation  given  in  Eq.  (2.6),  vie  compute  . 

veTys[<v;+vi)covf-<aiai+Vr),ln  I]  (4*64) 


This  expression  can  easily  be  evaluated  once  the  transmitter  and 
receiver  polarizations  and  the  bistatic  angle  are  known. 

The  covariance  of  the  channel  output  can  new  be  investigated. 
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Since  the  ec.tteeer.  are  i„  motion.  begin  by  exeminlng  the  Ho. 

«l  the  variables  x  =  (r.7>  in  the  phase  space  X.  Consider  two 

elements  dN(x.  t)  and  dN(x'.t').  Suppose  that  x  flows  into  x  as 

1 

time  Coe.  from  t  Co  T  and  flow,  into  *•  a.  time  goe.  from  t  to  V. 
where  t  =  {t+t’)/2.  For  sufficiently  .mall  values  of  t-t'.  we  may  put 

and 

**,  t'*tN  • 

-  “  xi+x2C-r)  w.m> 

where  x  i.  the  rate  of  change  of  the  position  vector  x  in  the  one 
ecatterer  phase  space.  Then 

w  -  a  fwft- «^^dN^.r, 


^**re  *  1  *1  (* j*  t )  and  Xj  -  *)•  It  follows  from  the  aes«ovip« 

tion  of  Poisson  statistics  tint 


Cov[w,w']  *  j^E  w  [x-*(^]w[;+  * ^  >  **  (4.471 

e 

where  we  have  assumed  that  the  velocity  field  x(x,t)  is  statistically 


i 

i 


*W5*Jy.*£W  ^mi  fr*.  it  ' 


J 


I 


so 


independent  at  the  measure  N(S,  t). 


Using  (4.^7)  and  making  assumption  (4.60),  we  obtain 

•  * 

Cov  [  K(t(w)(K*(t'(«')]  =  A*  exp  [  -i(»-«')to] 


"  A  *  exP  £  ~M®-«')to]  \  *XP  [  j  («-®')(aT+  aR).  7‘  ] 

E  (exp  [- (•  +  •')  ^~^.(5T+iR)*v(r’,T)] 
c[7'.v(7'ir)(^lG*[7'+v(7'J)^)l}<  n(7',7) >  d*7* 

jte  Q[*»y **[*'•* *  y(r!r)]£tr> dT  (4i68) 


where  v(r,t)  =  r(r, t)  and  we  have  assumed  tkt  v  and  y  are  indepcn* 
dent.  Let  F(v;  r,t)  be  the  probability  density  that  a  scatterer  has  a 
velocity  v  at  ( r,  t),  Eq.  (4. 68)  may  be  written  as 


■  ■  ■  m,  ■  will  pau. 


•I 


Cov  [ K(t, a), K*(t', »')]  *  A‘  exp[.i(«.«')L] 

o  w 

Jyeq [*•  *  -*  (rr  )]  Q* [ •* •  V  *  y  (rr )]  £(*> dt 

,J  j  {  exp  [  £  (»-«')(»T+  iR)  •  r*  ]  }  F(7;  7'.7) 

{g[  7'-7(7’.7)  (^)lc*  [  7+7(7 •  J)  (*-*)] 

exp  -  (»+»')  *r)‘7(7',7)1  |  <n(7',7)  >  d37d37* 


We  assume  that  Q(a,y)  is  sufficiently  wideband,  so  dial 


Q(w,y)  f«  Q{«  ,  y )  ever  the  signal  bandwidth,  where  s  is  the  angular 
©  © 


*  yj  _  _ 

carrier  frequency.  In  addition  we  assume  tint  —  (t'-t)B  <  <  1,  where 

© 


B  is  the  signal  bandwidth.  Substituting  (4.69)  in  (4.  97),  we  obtain 


Cov  [h(ti.h(f)] 


+• 


<2") 


W(nw.f-t) 


J  F*(»')  exp[-i*'  [ <'-to+  £(*t+  »r)[7’  +  7(7',7)^^~  *•'} 


g[7'-7(7'.7)  ^)]  c*[7‘*7(7',r»  (^)] 


F(»;i',t)  <  n(r',  t)  >  d3  v  d3  r’ 


.  lUU.illiJJ^  U9 


f 

& 


where 


% 

,  g -.  -  -v.?ff**v»,*»5s.-  m.rw^wv  --•  .?■ 
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w‘v*'->  *  Le  Q[v^K1ri)]  o'Cv^Ot)]  f<?> dT 


(4.71) 


Again  the  method  and  results  given  above  are  a  slight  generation  of 
those  given  by  E.  J.  Kelly  in  Reference  1 14]. 


For  the  dipole  scattering  example,  the  y -parameters  include 

* 

5»  TJ#  and  t  (see  Fig.  2),  where  the  angles  |  and  TJ  define  the  orien¬ 
tation  of  the  rotation  axis  and  (  -  a»^t+a  is  the  instantaneous  orienta¬ 
tion  angle  of  the  dipole  with  respect  to  the  u-axis.  Here  v^s  t  and  a 
is  the  initial  orientation  angle.  If  we  assume  that  every  orientation  ol 
the  rotation  axis  is  equally  probable  and  that  the  initial  orientation  of 
the  dipole  axis  is  uniformly  distributed  in  the  interval  0  <a  <2n,  we 
find  { 


where  t  =  t'-t  and  $  (r)  is  the  character!  Stic  function  of  p(v  )  defined 

•r  T 

ir.  Chapter  2. 


Sine  a 


and  f(t)  is  real  so  that  F*(-w)  =  F(»).  Equation  (4.70)  may  be  written 


u 


Cov  Ch(l),h(t*)]  *  A*  W(n^,t'-t| 
1 


J  i  f  (*-V  c  <V  V  [  r--v(7'.r,^)]} 

<V  V-  ['•♦*•»•.?>  0**)} 

F(v;  r',t)  <r  n(r',t)>  J*v  d*7‘ 


(4.751 


Now  f(t)  can  be  written  in  the  form  of  a  modulated  carrier, 

•  • 

f(t)  =  m(t)coi[»ot+T<>(t|]  (4.741 

Substitute  (4.  74)  in  (4.  73) -and  assume  that  m(t)  and  <p  (t)  vary  only 

slowly  compared  to  cos  n^t.  Then,  by  keeping  t'-t  sufficiently  small, 

we  may  drop  the  term  —  (?  +  a 

c  T  K 

of  m(t)  and  ^(l).  We  then  use  '.he  sum  and  difference  identity  on  the 
product  of  cosines  and  drop  the  sum  term  because  of  its  rapid  varia¬ 
tion  uith  r'.  Equation  (4.73)  becomes 


_  _  —  ,  t*-f\ 

)•  v(r,t)^-y j  from  the  arguments 


M 

Cov  th(t).h(f»  »  i  A*  W(n,  k'-t) 

CO  O' 

c*[  T'.v,?.t,(^)] 

mCl’lo+  £  <  V  V  • 7,1  mt  k’’ko+  £  <  V  V* 

cos  [  •~(iT+aR)-v(7'.tHt'-t)T  <  n(7\7)>  F{7;7',7)dv  dr* 

(4.75) 

The  variance  of  the  channel  output  can  be  deduced  directly 
from  {4.75).  Thus 


I 


i 

1 

j 


Varh(t)  "  Cov  [h(t),  h(t)] 

*  \  Al  ^™o^  c  ^*T+  *R)*^' ^  <  *&•*)>  d3' 

(4. 76) 

sphere 

w(u,o)  =  <  |G(wo.y)|2>  =  J  |Q(»0.y)|2f(7)dY  (4.77) 


The  quantity  j  m  (t)  is  proportional  to  the  average  power  being 
transmitted  at  t.  Likewise,  Var  h(i)  is  proportional  to  ine  average 
power  in  the  fluctuating  component  of  h(t),  received  at  t.  Recall  that 


X2  G_  G 
T  R 

..  ,3  2  2 

«")  rTrB 


o 


hence,  the  quantity 


•S 

w(.o>  jMg(7')|2  <  n(7,.t)>  ds7* 

has  the  significance  of  channel  cross  section. 

For  the  case  of  backscattering,  both  a  and  a  art  in  tka 

X  R 

direction  of  the  antenna.  However,  the  antenna  may  have  different 
polarization  for  transmitting  and  receiving.  If  one  introduces  a 

coordinate  system  (x,y,z)  with  the  z*axis  in  the  negative  (£  ♦  £  ) 

X  R 

direction,  and  assume  that  the  antenna  gain  is  nearly  uniform  over 
the  cloud  of  scattcrers,  then  the  cross  section  per  unit  range  is 
given  by 

o(z',t)  =  W(*o)  Jj*  <  n(x',y',*’;t)>  dx’  dy*  (4.7*1 

This  represents  simple  incoherent  scattering  from  a  collection  o t 
independent  scatterers,  i.e.,  the  total  cross  section  is  equal  to  the 
average  cross  section  per  scatterer  times  the  mean  number  of 
scatterers  illuminated  at  time  t* 

For  the  dipole  scattering  example,  we  have 

4 

8  —24*|ai(eo)i2<C2+U2+L2> 
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Chapter  5 

SCATTERING  BY  RANDOMLY  VARYING  CONTINUOUS 

MEDIA  WITH  APPLICATION  TO  COMMUNICATIONS 

' 

5. 1  Introduction 

Interest  in  the  problem  of  establishing  long  distance  commu¬ 
nications  via  forward  scatter  has  motivated  the  study  of  scattering 
from  a  perturbed  region  of  space  whose  electromagnetic  properties 
vary  randomly  with  spac£  and  time.  In  ti^  previous  chapter,  wo 
studied  the  properties  of  a  "discrete  scatterer"  channel,  in  which  the 
scattering  medium  is  made  up  of  a  collection  of  metallic  scattered* 
randomly  distributed  in  a  region  of  space.  The  scattering  charac¬ 
teristics  were  expressed  in  terms  of  the  electric  and  magnetic 
polarizabilities  of  each  scatterer. 

In  the  present  chapter,  the  properties  of  a  "continuum" 
channel  are  examined.  The  constitutive  parameters  of  the  medium 
are  assumed  to  vary  randomly  with  both  position  and  time.  The 
scattering  characteristics  are  expressed  in  terms  of  the  electric  and 
magnetic  susceptibilities  of  the  medium.  General  time  dependence 
is  assumed.  An  expression  for  the  channel  function  is  derived  in 
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term*  of  the  scattering  properties  of  the  medium .  The  approaches 
being  taken  here  are  parallel  to  those  of  Kelly  [  143 ;  howmr« 
slightly  different  techniques  are  employed  to  derive  the  scattered 
fields.  Also,  the  effect  of  antenna  gain  is  explicitly  taken  Into 
account  in  order  to  allow  for  the  fact  that  scatterers  may  flow  ia  and 
out  of  the  volume  illuminated  by  the  two  antenna  beams.  The  theory 
is  developed  !n  a  general  manner  so  that  it  may  be  applied  to  the 
study  of  a  wide  class  of  problems;  e.  g. ,  scattering  by  a  turbulent 
and  inhomogeneous  atmosphere,  ionised  gases,  and  ’  .bulent  wakes. 

5.2  Scattering  by  Continuous  Media 

Consider  ^  scattering  volume  which  is  filled  by  a  random 
continuous  medium.  The  interaction  of  matter  and  the  fields  ara 
described  by 

D  =  «  E  *■  P  <5. 1) 

9 

H  *  — -  B  -  M  (5.2) 

where  P  and  M  are  the  electric  and  magnetic  polarisation  vectors, 
respectively.  Outside  the  scattering  volume,  we  have  free  space, 
and  the  polarization  vectors  vanish.  The  electromagnetic  fields 
everywhere  are  governed  by  Maxwell's  equations  [S7]» 


The  aoluiton  can  be  written  as 


(9 


5(7.  t) 


where 


~  {Jj  P(7'.t')G(7,t;  7‘.t'>  dV  dt* 

o 

t* _ 

+  J  J  J  d*"t  V*KM(r,#t")lGCr,l;7,,lW,dl,J 

(5. 10) 


G(r.t;  r*.  t') 


Cs.ii> 


it  the  retarded  Green's  function  for  the  scalar  wave  equation  (Appen° 
dix  IV).  The  Herts  vector  for  the  primary  field  is  given  by 

V 

5o(r.t)  *  ~JJ[J  J(r'9t”)dt"lG(7.t;7'.t')dVMf  (5.12) 

o  "o 


Assume  that  all  sources  which  produce  the  primary  field  are  external 
to  the  scattering  volume.  Then  from  (5. 10).  the  Herts  vector  for  the 
scattered  field  is 


n#(7.t)  =  f"{JJ  P(7'.t')G(7.t;7'.t')dV'di' 

O 

te 

+  JJ[J  V'xM(7',t")dt"  jG(7,  t;  7',t')  dV'dt’^  (5.1J) 


where  the  domain  of  integration  is  the  scattering  volume  in  which  P 
and  M  differs  from  xero.  From  the  vector  identity  V'  X  04  » 

G  V'  x  M  -  M  xV  G,  we  have 
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f  (nxM)C  dS'  =  f  GV'xMdV  -  f  Mx?'G  dV  (5.14) 

JS  JV  JV 

where  the  integral  on  the  left  is  extended  over  the  surface  enclosing 
the  scattering  volume,  and  n  is  the  unit  outward  normal  to  this  sur¬ 
face.  Assume  that  the  scattering  volume  remains  within  a  finite 
closed  surface,  which  is  made  large  enough  so  that  the  surface  into* 
gral  in  (5. 14)  vanishes.  Equation  (5.  13)  can  be  rewritten  as 


V?.  0  *  ~  Jj  {  P(7'.  t')G(r,  t;  7*. f) 

+  j  M(7',t")xV'G(7,t;  r'.fjdt"]  dV’dt'  (5.15) 
Jo  J 

From  (5. 7)  and  the  relation  VG(  r,  t;  r',t')  =  -V'G(r,t;  r’,t'),  the 
scattered  field  becomes 


E#(r.t)  *  ~  jj  [P(7',t')*  V}VG(7,t;  7',t')  dV'  dt' 

'0 

*  Uo  fl  J  J  • t,)G  < t;  *'• V)  dV’  dV 

-U0—  JJIj  M(7,,t")xV'G(7,t;7',f)dt,'}dV,dt‘ 

dt  e 


(5.16) 


Let  the  origin  of  the  coordinates  be  located  within  V  and 
evaluate  the  scattered  field  only  in  the  far  sons,  then 
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|  r-r*  |  m  r  -  I  •  r* 


G  *  —  *(t'-t  ♦  -  -  -  *  •  7'| 

4nr  '  c  e  r  ' 

M  x  V'G  (a  xM)6’(t'-t  +  ---  5  .  7*1 

4ncr  r  c  e  f  1 

[P-Y'lV'G*  — ^[P-a?a  6"(t'-K4---  a  -*'l 
4«cZr  r  r  e  «  * 


where  is  a  unit  vector  along  the  redirection  and  6'  indicates 
differentiation  with  respect  to  its  argument.  Substituting  then* 
quantities  in  {5. 16)  and  making  use  6f  the  formula, 

j  f(x)  6(n,(x-a)  dx  =  {-if  ffcV> 

we  obtain 


E(r.t) 


.  2 
4nc  r 


+ 


1  • 

—  ax 

c  r 


where  P  = 


•  • 

P  and  M 


M  | 

X  j  [a  X  P(r',t- -  a  •  7'|T  dV* 

J  r  c  c  y 

w 

•  •  m 

j  M(r*.t-|*i£r-7*)dV’}  (5.1?) 

V 

a2  - 

*  — v  M.  The  scattered  field  can  be  inter- 
»•* 


preted  as  arising  from  a  distribution  of  induced  electric  and  magnetic 
dipole  moments  per  unit  volume  of  the  scattering  medium.  The  solu¬ 
tion  will  depend  on  the  manner  in  which  the  polarisation  depends  oo 
the  external  field.  If  we  assume  that  the  dependence  of  p  and  m  on 
E  and  H,  respectively,  is  linear,  we  can  write 


n 


P(r,l)  *  «o  J  X(7,t.»)‘  E(7,t-s)ds  (5.1i) 

M(r»l)  =  j  y  (7,1.*)  •  H(7,t-s)ds  (*.19) 

where  X  and  are  the  dyadic  electric  and  magnetic  susceptibilities 
of  the  medium.  Equation  (5.  18)  and  (5. 19)  also  express  the  fact  that 
the  induced  polarization  must  not  precede  in  time  the  field  which 

produces  it.  If  X(i,  t,s)  and  ^(r,t,  s)  vary  slowly  compared  to  the 

•  • 

primary  field  variations,  them 


P(r.t)  «  «o  J  £(7,t,*)»  E(r,t-s)ds 

© 


M(r,t)  «  J  y,  (r.t.s)  •  H(r,t~s)  ds 


Further  let  us  assume  that  the  presence  of  the  scattering  medium  is 
a  small  perturbation  of  free  space  conditions  so  thet  the  scattered 

A 

energy  is  small  compared  to  that  of  the  primary  field.  Then  we  can 
use  the  Born  approximation,  replacing  the  total  fields  E  and  H  by 
the  primary  fields  E^  and  Hq  within  the  integrals.  With  these 
assumptions,  Equation  (5. 17)  becomes 


F.ir.t)  * 

4«c  r 

»rx  J  |  +  ifr.7*-.ndVd.  (S. 20) 

2*  1  IB 

where  we  have  replaced  (t-  -  +  -  a^*  r')  by  t  in  the  argument*  of  X 
and  since  the  medium  is  slowly  varying.  This  is  the  same  as 
Kelly's  result  except  that  here  y  and  y  are  dyadics  instead  of 
scalars. 


K*U“, 


X  Cx<7'.  t. *)•  E  (7*, t -£♦- £  T'-stfJdYYls 


5.3  The  Channel  Function 

The  channel  function  for  the  continuum  case  can  be  derived  la 
the  same  manner  as  for  the  discrete  scatterer  case.  However* 
because  of  the  similarities  in  the  formulism  between  the  two  types 
of  scattering,  a  direct  analogy  can  be  established  so  that  the  results 
developed  for  the  discrete  scatterer  case  can  be  applied  immediately 
to  the  continuum  case.  Comparison  between  (4.2Z)  and  (5.20)  shores 
that  the  two  expressions  arc  equivalent  if  we  put 

a(x)  6{s)  dN(x .  t)  *— •  X(7',t.o)dV* 

“(7)  6(s)  dN(7,t)  —  jJ(7'.t.s)dV* 

Thus,  from  the  results  of  Section  4.3.  Chapter  4.  we  find  that  the 
channel  function  for  the  continuum  case  is  given  by 


*  ^*79***$  >**~. 


5 
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-i>H  +i*(a  ♦  £  )•»• 

K(t,«)*Ae  °f  C(T')Q(7M,a)e  C  T  R  dV 

°  V 


where 


,  G^-  G  X 

2  OT  OR 

•  *  2  r 

(4T,)  '*  'R 


Q{rt(t»)  =  ••"  v  [  e  •X(r,t,w)*  e  +  h  «  Y(r,t,®)*  h  ) 
VdifTc  R  T  R  T 


X(r.t.®)  *  J  X(r,t,s)  e"1"*  ds 


Y(r,t,«)  =  J  U>(r,t.«)  e"lw*  ds 


(5.21) 

(5.22) 

(5.23) 

(5.24) 

(5.25) 


Other  quantities  are  defined  in  Chapter  4. 

• 

The  condition  of  realizability  requires  that  the  polarization  be 
zero  for  s  <  0.  .  This  implies  that  X{  r,  t  ,«)  and  Y(r,  t,w)  must  be 
analytic  in  the  lower  half  of  the  complex  u-plane.  Since 
tQ-  -■  (a,y+  a^)  •  r'  >  0  for  any  point  in  the  scattering  volume,  it 
follows  that  K(tvu))  is  also  analytic  in  the  lower  half  e-plane. 


For  simplicity,  we  restrict  our  discussions  to  a  scattering 
medium  which  contains  only  a  single  constituent,  since  the  generali¬ 
zations  to  the  case  for  which  the  medium  contains  several  constituents 
are  straight  forward.  For  a  continuous  medium  with  one  constituent, 
having  number  density  n(r,t),  we  may  put 
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X(r,t,e)  =  X(n)  n(r,t) 


Y(r.t.«*)  =  Y(»)n(7 ,t) 


where  X(tu)  and  Y{«)  describe  the  dispersive  properties  at  the 
medium.  Thus,  the  channel  function  can  be  written  in  the  form 


-lot 

K(t,®}  =■  A  Q(®)e  ° 


i-(a  ♦£  U7* 

]  G(r')e  C  T  R  n(7',t)  dV  (5.26) 
* 


where 


Q(»)  =  X(®)»  e  +h  -  Y(n)  •  h  ) 

Y?Sc  R  T  R  T 


(5.27) 


As  an  example,  we  consider  the  scattering  by  an  ionised  gas 

(plasma).  We  assume  that  the  medium  is  non-magnetic  so  that 
as  — 

Y{r,t,es)  =  0.  The  electric  susceptibility  takes  the  form  [40} 


-dr  =&$ " 


(5.28) 


where  I  is  the  identity  dyadic,  e  and  m  are  the  electron  charge  end 
mass,  n(r,t)  is  the  electron  density,  and  v  is  the  collision  frequency. 
Thus, 


QW  =  >„4c2 


(*.*»» 


and  the  channel  function  becomes 


■  agflrtrJTSnsmr.’sss 


K(t, «) 


A  2  .i«t 

o  e  /  •  >.*  •  %  l"*o 

’  V4„  c2  m€o  *R‘  *t)* 

i-(Of_).  7* 

j  G(r')  e  C  n(r^,t)  dV 


(5. 30) 


The  channel  function  can  be  computed  if  the  density  profile 
is  known.  For  example,  let  us  consider  a  hypothetical  model  of  the 
plasma  for  which  we  take  to  be  a  spherical  column  with  a  density 


-a(t) 


t7*i* 


n(r’,t)  * 


J  a(t)  >  0 


(5.31) 


where  r*  is  a  vector  measured  from,  the  origin  within  the  plasma 
column,  a  is  a  constant  which  has  the  dimension  of  length,  and  a{t) 
is  a  time-varying  parameter'  which  defines  the  instantaneous  stse  os 
the  plasma  column.  For  simplicity  we  assume  that  the  antenna 
illumination  is  uniform,  G{?*)~  1,  and  that. the  plasma  is  collision* 
less,  v~0.  With  these  assumptions,  (5. 30). becomes 


A  •* 

A  • 

'W*  • 

JJJ  *  C  C  T  R  dx'  dy*  ds* 

(5.32) 


where  «*Z  =  n  e^/m*  .  In  (5.32)  we  have  assumed  that  ft(t)/a*  It 
op  o  o 

large  enough  so  that  n(r',t)  is  a  rapidly  decreasing  function  of  |  r'j 
and  the  integration  may  be  extended  from  - •  to  +  •.  Sine* 


r'  =  a  x'  +  £  y'  +  f  a* 
x  y  a 


and  a  and  a  an  unit  vectors  from  the  origin  to  the  transmitter  and 

I  A 

receiver,  respectively,  so  that  we  have 


a_  =  a  sin  6_  cos  0  +  a  sin  0_  sin  0_+  £  coe  0_ 
Tx  T  Ty  T  T  a  T 

i  =  a  sin  0_  cos  •  +  a  sin  0_  sin  #  ♦  a  coe  •_ 
Rx  R  R  y  R  R  a  It 


Equation  (5. 32)  may  be  rewritten  as 


A  V  -i«t_ 

K(t,w)  - - ?  (e_«  e_)  e 

V?n  CZ  R  T 

r+*  -  y,2+  *'2> l?  <V'*  V‘+  V'5 

iffe  *  c  y  dx' dy' da'  (5.3 


where 


=  (sin  0T  cos  #T+  sin  0R  cos  #R) 


=  (sin  0_  sin  0_+  sin  0M  sis 

T  T  R  2 


\  =  (cos  0T  ♦  coe  0Rt 


Carrying  out  the  integration  in  (5.33),  we  obtain 


■* 


F{»)  *  exp  C-«Z/B2] 


I*.  3*1 


Using  (5.36)  and  (5.37),  the  channel  output  becomes 

B 


h(t)  s- 


-  3  2 

A  /n  s  • 

Jc  (l  +4^coe2! 


171  1 V  V 


-2a(t). 


2> 


exp 


{-? 


B 


«  «M 


«-/} 


(S.3») 


Thus,  the  transmitted  pulse  is  widened  and  distorted  as  a  result  of 
the  electron  density  fluctuations  in  the  plasma  column. 

A  special  form  of  a(t)  may  be  represented  by' 


a(t)  *  a  ♦  a  b(t) 
o  1 

where  jb(t)|  <  1  and  a^b(t)  fluctuates  about  the  value  *lf  the 
time  variation  is  small,  (a 0^aj)  K<  i*  the  channel  output  is  approad* 


mately  given  by 

A  /n  a^B  • 


5.  * 


h(t}~ 


2 

4  yc  (t 


57^(V;t'[s-C^X> 


•*{- J^[*  * Ct X‘  -^>'“>1  «'  */}  <*•”» 


v  2  J 
9  Y 


•»  ^  JSW?»3fe»c*»-- 
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^•rc 


_  ?'_2  -j 

a  *  -g— -  cos  | 

c  a  * 

e 


in  the  stationary  case,  b(t)  =  0,  (5. 38)  reduces  i* 


b(t) 


A  /n  a^B  •  * 

o 


777737?  (v  V  “p  {-  i  ^  t'-*/}  cs.40} 

*  0  Y 


$•4  Statistical  Properties 

Some  of  the  statistical  properties  of  the  scatter  channel  are 
examined  in  this  section.  Again,  we  will  restrict  our  analysis  to 
that  of  a  single  constituent,  described  by  a  number  density  n(r,t}. 
We  assume  that  the  medium  is  randomly  varying  so  that  n(r,  t) 
becomes  a  random  process  in  space  and  time. 

From  Chapter  4,  we  recall  that  the  input- output  relation  for 
the  channel  can  be  written  as 


_  1  r 


M<)  =  ~  j  K(  t,  *)F (<b)  e1**  di 


where 


K(t,«)  =  f  k(t.s)  e”1**  ds 


(5.41) 


and 
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F(«)  *  j 

We  aiiune  that  f(t)  la  non- random  at'd  k(t,»)  la  real  a©  that 
X*(t,  -•).  The  mean  of  the  channel  output  la  gl  ran  hp 


E  h{t)  =  ^  j  E  K(t,e)  F(e)  dn 


(5.  42) 


From  (5.26),  we  have  Immediately 


E  K{t 


,.).A  Q(«).‘lo'°  j  0(»'l«  *  T  * 

O  •  u 


>  dV* 


C5.4JI 


where  <  n(r.t)  >  *  E  n{  r,  t)  la  the  mean  number  denalty. 

The  covariance  of  the  channel  output  la  given  1 

Cov  [h(t),h(t')3  *  —^2  J  J  Cov  [  K(t,«)»K*(t'»»')1 

(?n) 

exp  [■Hwt-ta't’JlFWFfa')  d®d»*  (5.44> 


Using  (5.26),  we  obtain 

Cov  [K{t.»).KV.»'»  =  A*Q<«)OV)  «P  C  -lU-rfH.] 
J  j  CU)C(7'>  exp  [  j  (»T+ 

Cov  ln(r,t),  n(r',t*)l  dV  dV 


(5.45) 


i02 


Tboi,  to  compute  the  covariance  of  h(t),  »e  need  the  covariance  ot 
n(r«t).  We  suppose  that  the  scattering  is  caused  by  inhomogeneities 
In  die  density  of  the  medium  in  random  flow.  In  order  to  facilitate 
die  computations,  we  shall  make  certain  simplifying  assumptions 
regarding  th*  '’ee.  Let  t  s  (tft')/2.  We  like  to  relate  n(r,  t)  to 
n(rj,t)  and  n(r',t')  to  njr^tt),  where  r^  is  the  point  to  which  7 
flows  as  time  goes  from  t  to  t,  and  r^  flows  into  r'  as  time  goes 
from  t  to  t'.  We  assume  that  t-t*  is  small  and  ignore  effects,  other 
than  the  flow,  which  can  change  the  density  of  the  constituent  being 
considered.  Then,  by  means  of  an  equation  of  continuity,  we  have 

n(7  ,7)dV.  *  n{7.t)dV 

(5. 46) 

n{72.7)dV2  *  n(?*,t')4V* 

where  dVj  and  dV^  are  two  elementary  volumes  at  r^  and  r^,  respec¬ 
tively.  We  also  assume  that  the  velocity  of  the  random  flow  is 
independent  of  the  density  n{  r,  t).  Then,  for  small  values  of  t-t*, 
we  may  writo 

V  '+v<vo(£ji) 

<5.«> 

-  r.  r  t‘-t  \ 

where  v(r,t)  is  the  velocity  of  the  flow  at  (r,  t).  Thus,  if  #{r)  Is 
some  function  of  position  and- if 


<P  s  f  *(7)n(7,t)dV 

“V 

x  J  #'(T')  n(7',l')dvt 
v 

^,ea»  (rom  (5.46}  and  (5.47),  9  and  <p'  can  also  be  expressed  an 

*  *  L  '[v'iv‘|Ct^»iv'1'1vi 

vl 

’  =  J„  *’[  V  7(72'7>  C'Ti)]"‘72-*'>  dV2 

V1 

It  follows  that 


Cov  [  9t <p'  1  *  f  j  t(r)6'{r')  Cov  [  n(  r,  t),n(7'#  t’)}  dV  dV' 

TV 

1  jv  1/  ,[7i-7‘7rr»Cir1)]»,[72+7«72*r»Ciyi)] 


1  ”1 


Cov  ln(7l,7),n{7l§7)')dVldVl 


(5. 48) 


In  this  manner,  we  may  put 


-  -  -  -  r- 
Cov  [n(r1#r),n(72#7)]  =  aQ — ^ - ,  VVO 


C5.49J 


Applying  (5.48)  to  (5.45),  we  obtain 


8HP 
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Cov  r.K{l,®).K*(lJ»')]  =  A*Q(*)Q*{*‘)  exp  C-i{«-*')|  ] 

L  L  ,xp  [*  c  (V  *R>- ”t>] 

1  2 

Ec[7rv(7I.r)^)l  G*[V7{V7,(rr)] 

«xp  {- l-  «»x+  *R>  •  C«'7{72 .7)  +  •7(7I,oi(l~i)J 

A  C  "V^  •  V  Vr)  dvi  dV2  (5*  50) 


If  «•  della* 


Vr2 


and  p  *  *2-  7j 


Equation  (5.  50)  may  be  written  a«  ’ 

Cov[K(t.»).K*(f.«')?  =  A^Qtui^^JexpC-Ka-a'Hj 

Jj  <sp  {  ^  (*t+*r>*  [<•-•')*  -(•+•') f  1 

E c[7 - f  -  y{7- §  .rX^)]o*[74 *7(7, | .T)(^)] 

"p  {-  H^rDCrr)' V  *‘r>-  [  *«*'♦  I  .o- v*t-|  ,r>  1} 

«<*  {-  7  C-t^Ct5)  <  v  *•„>•[  ** !  ••>  *  «<-!  .«>] } 


A(r,p»t)  d*  r  d*  p 


(5.  51) 


We  assume  that  A(r,p,t)  varies  slowly  with  r  and  vanishes  relatively. 

rapidly  as  |  p  J  increases.  Specifically,  we  assume  that  A(7,7»7) 

becomes  negligibly  small  when  (  p  (  exceeds  some  correlation  length 

p  which  is  small  compared  to  the  modulation  wavelength  X  *  2nc/B» 
c  m 

where  B  is  the  signal  bandwidth.  This  implies  that  the  signal  illum¬ 
inates  many  uncorrelated  scattering  volumes  with  nearly  uniform 
modulation.  Thus,  by  keeping  p  sufficiently  small,  we  may  assuma 


7(7  ±  ~  ,7)  «  7(7,7) 


Equaj^n  (5.  51)  becomes 


Cov  [  K(t.  w),  K*(t',  u*')1  *  A*Q(»)qV  )  exp  [-t(»-n')tj 

P  o 

G  [7.  |  -  7<7,7)  ; V [?  ♦  i  ♦  7(7,7, 

exP  (»  +  ■')  C^T"^*t+  *  7(7,7)  ^ 


F(v;  r,*t)  A(r,  p,t)  d*  v  d*  r  d*  p 


(9,  92) 


where  F(v;  r,  t)  is  the  probability  density  for  the  velocity  of  Idas  flow 

at  (7,t). 


Now  we  assume  that  Q(u>)  is  wideband,  i.  e. ,  Q|«)mQ(«  ), 
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that  —  (t’-t)B  <<  1 ,  where  ^  is  the  angular  carrier  frequency  and  B 

*  • 

is  the  signal  bandwidth.  Substitute  (5.  52)  into  (5.44),  we  obtain 

.  I 


cov  [hio.Mfn 


I  «<•„>  I 


JJJ  '{‘•vi<W-[7-|-’',7*r,Cir)]} 

‘  {*  -‘.+  7  <  V  *.»*[' 7  +  I  * 7l7-7>  (tT  )]} 
o  [  7  - 1  -  7(7.7)  (££)]  G*[7t|<  7(7.7)  (i^)] 


F(vi7.t)  A(r.p.t)  dN  d*  r  d5p 


(5.  53) 


where 


f(t)  *  ~  J  F{ml  e1"4  d» 


As  in  the  discrete  scatterer  case*  we  assume  that  f(t)  is  of 
the  form* 


5(t)  =  m(t)  cos  (  »  »  +  9_(t)3 

o  o 


(5.54) 


where  m(t)  and  9  (t)  are  slowly  varying  functions  of  time.  Substitute 
o 

(5.  54)  in  (5.  53)  and  assume  that  p  and  t'-t  are  sufficiently  small*  so 
that  we  may  drop  the  quantities 


1  7(aVV-[f  +  7,7-r,(iT)] 


-  xt'-t^ 


i 


.  ■ 


from  the  argument*  of  m(t)  and  9  (t)«  We  then  use  the  tom  and 

o 

difference  identity  on  the  product  of  cosines  and  drop  the  sum  term 
because  of  its  rapid  variations  with  r.  Equation  (S.  53)  becomes 


Cov  [  h(t),h(t')3  *  -j-  |Q(*o)|2  JJJ  F(7;7.7)A(7.?J) 

G  [  t  -  |  -  v(7. t )  C*[7t|  +  . O (^-)] 

m[  t-t#*  |  (»T+  »R)  •  r  }m[  j  {»,♦  fR)  •  rj 

CO*  1  "c^  P  +  v(r,t)(t'-t)3  j  d  vd  rd  p  (5.55) 


The  antenna  gain  patterns  can  influence  the  covariance  functiou 
greatly  when  the  value  of  v(  r#  t  )(t'-t)- exceeds  the  correlation  length 
t>c«  As  v(r,t)(t'-t)  increases  to  a  large  enough  value,  the  Integra* 
tion  will  be  carried  over  disjoint  volumes  of  the  scattering  materials, 
and  the  covariance  function  will  drop  rapidly  to  aero. 

The  variance  of  the  channel  output  can  be  obtained  directly 


from  (5.  55).  Thus 


10ft 

Varh(t)  *  Cor[ hfl),**®)) 

A*  2  -  - 

*  ~T  Jj  A(7.7.t)C(7.|)G*(T»|» 

m2[t-»oti<iTt  ;ri.T]cm  [-f-<»T+  »g)-Vj<i57fi5(r 

«*.») 

Since  jm2(t)  is  proportional  to  the  average  power  being  trans* 
mitted  at  t,  hence  Var  b(t)  is  proportional  to  the  average  power  la 
the  fluctuating  component  of  h(t),  received  at  t.  Substitute  the  value 
of  A*  from  (5.22)  in  (5.  56)  and  compare  with  the  bistatic  radar 
equation*  we  see  that  the  c  hannel  cross  section  is  given  by 

o(t)  =  |Q(«o)l2  JJ  A(7J,i)G<7-|)G*(7+f  ) 

cos  [  ^  <ST+  aR).  7]  dJ7d37  <5.  57) 

In  many  practical  situations,  it  is  possible  to  put 
A(r,p,t)  =  C(p)  Var  n(r,t) 
where  C(0)  =  1.  Then  (5„  57)  may  be  written  as 

o(t)  =  |Q(«0)j2  JjC(r  -  | )  G*(r  +  ^  )  Var  n(7,t) 

C(p1  cos  (aT+  aR)<  p]d37  d3  p 
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where  0  is  the  bietafcic  sagie, 


Appendix  I 

APPLICATION  OF  THE  LOR  ENT Z  RECIPROCITY 
THEOREM  TO  SCATTERING  PROBLEMS 


Let  (Ej.Hj)  and  (E^.H^)  be  two  different  distribution*  of 
electromagnetic  fields  of  the  same  frequency,  which  can  be  estab¬ 
lished  in  a  region  of  space  containing  no  sources.  Each  set  of  field 
vectors  satisfies  the  Maxwell's  equations.  Lorentx  reciprocity 
theorem  sta»'  . 


n*  (E,X  H_-  E.X  H.)dS  »  0 


a.i> 


where  n  is  the  unit  normal  vector  to  the  surface  S  enclosing  -’he  region 
of  interest.  Although  there  are  many  forms  of  the  reciprocity 
theorem  [32}-  [35],  the  form  we  use  here  is  that  due  to  Lorents 
[13],  [36], 

In  general,  the  scattering  geometry  consists  of  a  transmitter 
which  illuminates  the  target  and  a  receiver  at  which  we  observe  the 
scattered  signal.  A  schematic  representation  is  shown  in  Fig.  3.  In 
particular,  we  seek  an  expression  for  the  voltage  induced  at  the 
receiver  terminals  by  the  field  scattered  by  the  target* 


m 


Figure  8.  Schematic  Representation  of  Scattering 
Geometry 


IM**!!*?"'*. 


K*'«PS**5*i* S'.'  'C— 


Lei  (EjiHj)  be  the  fields  of  the  receiving  antenna  with  driving 

current  at  its  feed  point  with  no  other  sources  present;  i»e. ,  both 

the  transmitter  source  and  target  are  removed.  Let  (E^H^)  be  the 

total  fields  obtained  with  both  the  transmitter  source  and  the  target 

in  place.  It  is  assumed  that  the  same  current  L  is  delivered  to  the 

R 

feed  point  of  the  receiving  antenna  either  with  the  transmitter  source 
and  the  target  removed  or  in  place.  Thus,  E^  and  represent  the  ' 
total  fields  which  would  be  obtained  in  space  if  the  receiving  antenna 
were  used  as  a  transmitting  antenna,  while  and  represent  die 
sum  of  the  transmitter,  receiver,  and  scattered  field*. 


The  surface  integral  in  (!.  1)  can  be  written  a* 


j  »  •  ( E  X  H_  -  E.  xH.)dS  =*  8 

JS  +  S_+  S  ♦  S_  1  *  1 

•  T  •  R 


a.  2) 


The  integration  over  the  spherical  surface  vanishes  on  account  at 
the  radiation  condition  as  the  sphere  recedes  to  infinity  [3* ].  For 
simplicity,  we  replace  the  transmitter  by  an  equivalent  source  J^,, 
which  creates  the  incident  field.  From  the  divergence  theorem,  vre 

have 

J  n-CEj*  H2-  E^x  Ht)dS  =  -J  V  •  (EjX  H?-  E^x  Hj)dV  $1.3) 


Within  the  volume  V^,  Maxwell's  equations  are 


Ill 


7  X  E.  *  -i»M  H,  7  X  E  *  -ispt  H„ 

i  o  i  <  ex 

7xH*  H»tE,  7xH  »  +i«r  £„♦  J. 

1  o  1  2  o  2  T 

Eq.  (1. 3)  become* 

j  S*(ElXH2-  ^XHjJdS  *  J  E l*JT  (1.4) 

ST  VT 

The  integral  over  the  surface  of  the  seatterer  can  be  evaluated 
in  a  similar  fashion.  That  is,  we  replace  the  seatterer  by  an  equiva¬ 
lent  distribution  of  current  induced  by  the  incident  field.  Then, 
within  the  volume  V>f  enclosed  by  the  surface  S>#  the  fields  satisfy 
Maxwell's  equations 

7  x  E j  *  -isuJHj  7  x  E^  *  -isU^Hj 

7  x  H ,  ~  ♦  is  «  H ,  7  X  H,  *  +i*t  E,4J 

1  o  1  .2  o  2  s 

Again,  applying  the  divergence  theorem,  we  obtain 

j  n  « (EjX  H^-  EjXHjJdS  *  J  E^J^V  £.  55 

®s  V. 

The  surface  includes  the  exterior  shielding  surfaces 
{assumed  to  be  perfectly  conducting)  of  the  receiving  system,  tit® 
inner  metal  surfaces  of  the  antenna  and  transmission  line  down  to  t*® 
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feed  point  (shown  by  the  plane  AA*  in  F»g.  8),  and  the  cross- 

sectional  plane  AA*  inside  the  transmission  line.  Since  n  •  E.  x  H_  » 

12 

n  x  Ej®H^  *  0  and  o*EjX  Hj=  n  x  E^*Hj  =  0  on  the  metal  surfaces, ' 
die  only  n  on-vanishing  contribution  to  the  surface  integral  over  S  Is 
that  over  the  cross-sectional  plane  AA',  Thus, 


J  n.(EJxH2-E2xHi)dS  =  J  n°  (E^x  H^-E^XHjjdS  (I.$> 

*R  AA' 

The  evaluation  of  the  surface  integral  over  AA*  requires  some  const* 
deration  of  the  field  distribution  inside  the  transmission  line^  For 
the  present  purpose  we  consider  a  coaxial  line  [  13],  although  such  a 
choice  is  arbitrary  [36].  For  a  coaxial  l!ne,  the  fields  in  cylindrical 
coordinates  are  given  by 


E 


* 


r  ~  TfnfbTbT)  *  n$ 

o  t 


Ivr 


a.  ?) 


where  h.  and  b  are  the  inner  and  outer  radii  of  the  line,  and  V  and  1 
*  e 

are  the  line  voltage  and  current.  Substitution  of  {1.7}  into  (1.6)  gives 


'AA*  1  2 


SjSHjtdS  « 


V  £  -V  1 

ii  n 


a. ») 


Let  V  fee  the  total  voltage  induced  at  the  terminals  of  the  receiver  by 

the  transmitter  field  and  the  scattered  field,  then  V  *  V.*  V, 

2  1 


V  - 


MftdSt  «. 


ii  !>*  •>*  ^  - 


Appendix  Q 

PROBABILITY  DENSITY  FUNCTION  FOR  THE  RETURNED 
WAVEFORM  FROM  A  RANDOM  COLLECTION  OF 
ROTATING  DIPOLE  ?*,ATTERERft 


The  probability  density  function  for  the  echo  waveform  can  be 
evaluated  by  first  finding  the  characteristic  i metion  of  the  distribu¬ 
tion  [  24 1,  [251.  For  simplicity  we  let 


"» 

*  J  •  •  •  J  P<V*»0  y  •  -j  PCW.  P‘5ff^«r 

••  0  -• 


pfw)d«  *  ||  p(w  )dw. 

i*I 


«•  «4  «4 


whore  w  represents  a  ,6,0,. » .  etc.  Hence,  we  can  write. 


p(*s  * £  r. ..r  pta 

n*C 4  V  0 
6 


4T 

...dw  f. . .  f  p(  t)dt  8(x-x  )  (II.  1) 
:  v  r  j  4  a 

-T 


From  Eqt.  (Jt2)  and  (3,4|,  the  echo  signal  can  be  written  in  the 
form,* 
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! 


ntt'rWS 


1  IT 


x<*>s  £  Mi. 


M  #^‘*VLC  ««  V  0—  <W  ♦  L  .«  <#,.  V]  0,2) 


where 


V'-V  *  -r°<V  V*«k,<*-‘k» 

*0 

*1  *  <V  V"-‘k»  *  V-V*  V“V*k.»+x(«k-V 

*. =  V-r|,!4k,w,|k| 

The  characteristic  (unction  of  the  distribution  becomes 


♦<«)  *  j  I»fcJ 


•  dx 


«*»  - 

"  r  r  -  „  -  +T  n  vfti  ‘I 

*  J"V“# . •r^O-<**rj-J  TT  'T  d 

I  -x  ^*1 

a 

CXP  lU  kn*k<{#lk)t-C  COa  VUco*  <VV* 

-  {  v(t]M8 

„«0  sr{  j  v<f)dfJ...J  P<V--*r  1^0- • 

W| 

•  exp  [iu  cos  $t+  u  cos  (^4#^  4  L  cos  (0^*1  J)]  J* 

1  “p{  pl*o . 

wi 

exp  [iu  #k(t,f)(C  cos  $x+  u  cos  {^4  0#)  ♦  L  cos  (*j-##W]  } 


*■'**, -  ....  *.  -. . , — _  ^  «- > «■'  5r^  5 


lift 


Expanding  Iht  exponential  inaidr  the  integral  in  a  power  series*  w« 

obtain 


t(u)  *«xp{-f  v(t)dt  ♦  E  Erf  vU'Jdt' f« •  .f  p{aA,..„,»  5 
-T  l*0m*0  n*0  -T  0  r 

*  -r  lyin  a 

d*o* **d®r[ ,u  V‘.^]  Ti^f  cosm(#i+  #,> 

ft  +T  mm 

coon(ft  -#  )  *  exp  {  -  ~  o2  ♦  f  vU'Jdt'  EEC 

•T  1*1  m*l  n*I 


IM 


#• «m+n 


llmlal 


W(l,m,n;  t,t'|| 


%  2 

where*  for  simplicity*  we  have  set  <  x  >  s  0,  <  x  >  *  o  *  and 


W(l*m*n;t«t')  =  J**»J  p(aQ... . ,oBr)da0...d*r  £tk(t,t')l 


1+m+n 


•  Cl^nLnco.Slc«s,B|»1+  0#)cosn(#j-  0#) 


More  conveniently*  we  can  write  the  triple  summation  as 


*  *>•*> 

q*e  1  n 


where  q  *  1  +m*n 


Thus* 


tiu)  =  exp 


,  2  ,  +T 

r  u  _2 .  f 


[-Vo**  [  v(f>dt’  E  W  (t.f)) 

'  2  -T  q*3  qI  *» 


(E.J) 


Ilf 


By  ttking  the  lover..  tr.n.form  of  the  eher.cteri.tlc  fooctloo. 


obtain 


P(»l  ■  ^  J  ♦<«>  do 


•  h  J/--0 -  r1*]'  -p[ 0..4. 


«her« 


s  J  v(l')W (t,f)df 


01.5) 


If  we  expand  the  second  exponential  of  (II.*)  |n  ,  power  series,  we 


have 


*  , 
»  2 


i  r  Me  t  * 

■  si/- '•  «... 


Equation  01.6)  give*  the  time-varying  probability  distribution 
of  .he  echo  waveform.  We  will  show  that  the  probability  density  at 
a  particular  instant  of  time  tends  to  become  a  Gaussian  probability 
density  as  the  echo  rate  becomes  large.'  If  we  let 


°  *  0  and  V  *  «« 


Equation  01*  6)  can  be  written  as 
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ppj .  -L^  |Vi»“-  YheCillfi  t  (ly)4  °«  . 

2no*  2-  *-  3!  3  *  41 

O  0 


} 


w«  not*  that 


j  <lY)ViY°“  2  4y 


_  .n 

QQ  -  •» 


Thu».  ev.lu.tlng  the  integrals  .nd  putting  p(x,  ,  op(a).  ».  obl.[, 

2 

r 

•  «• 

P<*> 


_L_  .  *o* 

V2ii  o 


cn.7> 


wher«  Hk»re  Hermite  polynomials  defined  as  [3»3, 


A  «2 

SL  k  SL 

HkW  *  (-l)k  •  2-i-  * 

da* 


We  assume  here  that  all  moments  of  the  distribution  exist  and  note 


Ihst  a  can  be  written  as 


-  [C 


k 

v<t’)W2it,f>dt'J 


01.8) 


From  Eqs‘  tH-5)  and  01.8),  we  have 
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DYADIC  GREEN'S  FUNCTION  FOR 
THE  VECTOR  HELMHOLTZ  EQUATION 

The  vector  Helmholtz  equation  for  the  electric  field  with 
c/*"1  time  dependence  is 

V  x  V  x  E(7)  -  k2E(  r>  =  -lep  J(7)  (III.  1J/ 

2  2 

where  k  =  w  tie .  A  solution  for  th«  inhomogeneous  vector  wave 
equation  can  be  written  as  [28]# 

E(7)  *  -i unx  f  r(7,7').  j(7')dv‘  (iu.2) 

Jv 

where  the  double  bar  indicates  a  dyadic  quantity.  Substitution  of 
(ill.  2)  into  (III.  1)  gives 

f  C  vx  vx  r(7,7')-k2T(7s7')]  •  j(7')dv  =  7(7)  tii,3) 

Since 

j(7)  =  f  r*j;7’) 6(7-7') dv* 

where  I  is  the  identity  dyadic.  Equation  (III. 3)  becomes 
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f  C  vxvx  7(7,7'}-  k  7(7,7*) -76(7-7  )]*j(7)dv’  *  o 

Jy 

since  this  relation  holds  for  any  arbitrary  -U  r')e  r(r,r')  must  satisfy 
the  differential  equation, 

v  x  v  x  7(7,7)  -  k2  r{7,7')  =  "6(7-7’)  (ui.4) 


We  require  that  the  solution  satisfies  the  radiation  condition  for  r*<«*. 
The  r-dependence  must  be  in  the  form  of  an  outgoing  spherical  wave. 

•  —  —  2  — 

Using  the  vector  identity  V  X  VxA  =  V(V  *A)  -  V  A,  011.4)  may 
be  written  as 

(v2+  k2)  7(7,7')  =  v[v  •7(7,7')]  -  76(7-7')  on.  5) 


Taking. the  divergence  of  011.4),  we  find 


▼  •  r(r,7')  = 


Equation  0SI.  5)  becomes 


-  \  v  6(7-7’) 

k 


(V2+ 


) 7(7,7')  =  -^F+  ~  vv)  •«*-»') 


Let 


7(7,7')  =  ^F+  ~  vvj  g(7,7') 


011.  i) 


0U.7) 


US 


where  P  and  Q  are  two  vector  functions  of  position  and  n  Is  th« 
outward  unit  normal  vector.  In  order  to  apply  the  vector  Green's 
theorem,  we  scalarly  multiply  the  dyadic  equation  (III.  4)  cm  the  right 
by  an  arbitrary  constant  vector  a  and  obtain  a  vector  relation  f ©? 
r(r.r').  a, 

vxvx  r(7,7'j.  I  -  k2r{7,7ij.r=  16(7-7*)  (m.is) 

Let  P  =  T(  r ,  r')  •  a  and  Q  =  E(  r').  Equation  (III.  12)  together  with  Eqs. 
(HI.  1)  and  (UI.  13)  give 

E(r)*a  =  -iu>u  f  [  r(r,r')*»l*  J (r'JdV' 

JV 

-  j  n*.  {cT(7.7.).73xv*xE(7.)-E(7*)xv.xCn7,7.).r]}d3. 

s 

(in.  14) 

From  the  .Maxwell's  equation  V'  xE(r')  =  -imiJHfr1)  and  the  vector 
identities 

n'  •  [  7(7,7') -I]  x[V'  xE(7')3  =  -n'  x  [  7'  xE(7')].  [  r(7,7')»a] 
n' .  E(7')xVx[r(7,7!)*al  =  n'  xE(7')- V'  x  [ 7(r,r'>*a3 


Equation  (IU.  14)  may  be  written  as 


flur  .ww  jim  m*  ■ 
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E(7)*a  *  -ia»u  J  f.  r(7,7').I]*  J(7')dV' 


liwutfi-xHt^n*  Cr(r.r')*a]-[5IxE(7')].V,x[r(7t7*)-7lJdS» 


cm.  is} 

Since  the  constant  vector  a  may  be  chosen  arbitrarily,  Eq.  (III.  15$ 

(Ives 


E  (r)  =  -imp  j  7(7,7')  .  J(7')  dV 

V 

-  j  {imuCfl*  xH(7')].r(7,7')-[n' xe<7»)].  vx 7(7,7’)  } ds» 

8  (HI.  16) 


The  volume  integral  is  a  solution  to  the  inhomogeneous  vector 

«— ■  _ 

Helmholtz  equation  with  a  spatially  distributed  source  J(r').  The 

surface  integrals  are  solutions  to  the  homogeneous  Helmholtz  eq  uativu 
-  2  — 

VxVxE-kE  =  0in  terms  of  the  boundary  values  on  S.  The  first 
term  in  the  surface  integral  may  bo  interpreted  as  a  contribution 
from  a  surface  current  of  density  n  x  H  on  S,  and  the  Second  term 
may  be  considered  as  a  contribution  from  a  magnet;**  current  sheet 
of  density  E  x  n.  For  a  perfect  conduct  r-  the  tangential  component 
of  E  vanishes  on  S,  then  nxE  *  0, 
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GREEN'S  FUNCTION  FOR  THE  TIME-DEPENDENT 
WAVE  EQUATION 


Consider  the  inhomogeneous  scalar  wave  equation, 

vz*  ■  H: s  ‘f('*k)  (IV*l) 

c  dt 

where  f  ( r,  t)  describes  the  distrioution  of  sources  which  are  func¬ 
tions  of  both  position  and  time.  The  solution  [39]  -  [40]  may  bo 
written  as 

Hr.t)  =  f  f  G(r,t;7'.f)f(7\f)dSr'  dt'  (IV.2J 

* V4o 


where  the  Green's  function  G(r,t;  r'.t')  satisfies  the  differential 


€t|\Sntiv'n| 


V2c-  ~  =  -  6(r-f)  6(t-t') 

c  dt 


(IV.  3) 


To  find  G,  we  consider  the  Fourier  representations. 


G(r,  t;r’,t’)  =  Jd3k  J  dx  g(k,u>)  e 


-ik.(r-r')  «i»( t-t') 


(IV.  4) 


and 


6(r-?)  6(t-t’|  =  fd3k^dBe'k'(^■^,)e+il^,<t"t,,  (IV.  5> 

<2")'  4 
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*>™****^5m!mm8 
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Substituting  (IV,  4\  and  (IV.  5)  into  (IV.  3),  we  find 


g(k.J>) 


(2"> 


Ck*-4) 


(IV.  6) 


Hence 


C(7,t;T'..')x  -i-3  fd’kf - is- 

(2")  J  J  rJt  •  n 

ik'j)  (tv.v 


The  Green's  function  satisfying  (IV.  3)  represents  the  wave 
disturbance  caused  by  an  impulse  source  at  t=t\  located  at  r*r\ 
Thus,  G  gives  the  description  of  the  effect  of  this  impulse  as  it 
propagates  away  from  r=r'  in  the  course  of  time.  It  ca  t  then  be 
assumed  that  G  and  dG/dt  should  be  zero  for  t  <  t';  i,  e. ,  no  effect 
should  precede  the  ca.ise  in  time.  Furthermore,  for  t>  t',  the 
wave  disturbance  propagates  outwards  as  a  spherically  diverging 
wave  with  a  velocity  c. 

Equation  (IV.  7)  can  be  evaluated  by  application  of  the  Caucny 
residue  theorem.  However,  in  order  to  make  G  vanish  for  t<  t', 
we  imagine  that  the  poles  at  uu  =  +  ck  are  displaced  above  the  real 

v 

axis  by  an  arbitrarily  small  amount  c  in  the  complex  a-plane.  The 
integral  over  the  lower  half  plane  (t  <  t* )  will  then  vanish,  while 
the  integral  over  the  upper  half  plane  will  give  a  non -vanishing 


) 


l 

? 
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where  R  =  |r-r’  J  and  T  =  t-t'.  The  second  term  in  (IV.  9)  does  not 
contribute  because  of  the  conditions  of  realizability  mentioned 
earlier.  Therefore* 


G(r,  t;  t* ) 


4n  |  v-sr'l 


^t-t' 


t-t'  >0 


(IV.  10) 
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